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Foreward 


The  research  and  development  work  "Investigation  of  ice  dynamics 
in  the  marginal  ice  zone”  was  commenced  by  Dr.  Matti  Lepparanta 
on  January  24,  1983,  at  the  U.S.  Army  Cold  Regions  Research  and 
Engineering  Laboratory,  Hanover,  N.H.  The  work  has  been  done  in 
cooperation  with  Dr.  William  D.  Hibler  III  from  USACRREL  and 
consequently  this  report  is  co-authored  by  Dr.  Hibler. 

This  work  consists  of  a  theoretical  and  an  experimental  part. 

The  former  involves  dynamic  modeling  of  the  ice  flow  in  the 
marginal  sea  ice  zone  and  the  latter  ice  kinematics  measurements 
in  MIZEX-83  experiment  in  the  Greenland  Sea. 
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1.  Introduction 


The  marginal  ice  zone  (MIZ)  is  defined  as  that  part  of 
the  sea  ice  cover  which  is  close  enough  to  the  open  ocean  to 
be  affected  by  its  presence  (Wadhams,  1980).  It  is  a  region 
of  complex  atmosphere- ice-ocean  interaction  extending  to  100-200 
km  from  the  ice  edge.  Well-defined  MIZ's  are  typically  found 
in  the  European  sector  of  the  Arctic  ice  pack  (Fig.  1)  and 
along  the  northern  edge  of  the  Antarctic  ice  cover.  The  location 
of  the  ice  margin  has  a  significant  role  in  the  world  climate 
system  and  hence  the  MIZ  has  recently  become  a  region  of  intensive 
research. 

International  scientific  program  called  Marginal  Ice  Zone 
Experiment  (MIZEX)  was  established  in  1981  (Wadhams  et  al . , 

1981).  Its  objective  is  to  understand,  in  the  mesoscale,  the 
physical  processes  by  which  ocean,  ice  and  atmosphere  interact 
in  the  region  of  the  ice  edge.  The  main  field  effort  of  the 
MIZEX  program  takes  place  in  the  Greenland  Sea  because  most 
of  the  heat  and  water  exchange  between  the  Arctic  Ocean  and 
the  rest  of  the  world  ocean  occurs  through  Fram  Strait,  the 
channel  between  Greenland  and  Svalbard  (Johannessen  et  al . , 

1983a).  A  summer  experiment  will  be  carried  through  in  1984 
and  it  was  preceeded  by  a  pilot  study  in  1983. 

A  sea  ice  cover  modifies  the  air-sea  energy  exchange, 
penetration  of  solar  radiation  into  the  sea  and  the  salt  budget 
of  the  upper  ocean.  Ice  dynamics  deals  directly  with  the  kinetic 
energy  exchange  but  also  affects  the  thermodynamics  indirectly 
through  opening  and  closing  of  leads  and  advection.  In  turn, 
freezing  and  melting  change  the  strength  of  the  ice  which  determines 
its  resistance  to  differential  motion.  Thus  dynamics  and  thermodynami 
are  coupled.  However,  on  short  time-scales  sea  ice  drift  can 
be  treated  as  a  purely  dynamical  problem. 

The  viscous-plastic  sea  ice  model  of  Hibler  (1979)  can 
simulate  the  large  scale  drift  of  the  Arctic  ice  pack  quite 


well  (Hibler,  1979,  1980b;  Hibler  and  Walsh,  1982)  but  the  model 
resolution  has  been  too  cosrse  for  the  characteristic  scales 
in  the  marginal  ice  zone.  The  model  has  been  applied  to  the 
East  Greenland  area  on  a  finer  (40  km)  scale  (Tucker,  1982) 
and  to  sea  ice  in  the  Bay  of  Bothnia  on  a  finer  scale  (18.5  km) 
yet  (Hibler  et  al.,  1983).  Although  these  scales  are  still 
a  bit  coarse  for  the  MIZ,  the  results  indicate  that  the  model 
works  well  when  its  scale  goes  down.  Ice  velocities  near  the 
ice  edge  were  not  well  simulated  which  was  considered  to  be 
at  least  partly  due  to  an  improper  wind  field  and  improper  current 
fields.  However,  in  the  MIZ  the  ice  is  more  fragmented  than 
in  the  central  pack,  with  substantial  variations  in  compactness, 
and  the  validity  of  the  particular  viscous-plastic  rheology 
needs  to  be  examined. 

In  this  work  the  dynamics  of  the  pack  ice  in  the  MIZ  are 
considered  both  from  theoretical  and  experimental  point  of  view. 
Emphasis  is  laid  on  the  kinematic  properties  of  ice  and  on  the 
role  of  ice  interaction.  Modeling  studies  are  carried  through 
to  examine  how  various  standard-type  pack  ice  rheologies  may 
modify  the  MIZ  flow.  Field  data  on  ice  kinematics  are  presented 
and  analyzed.  They  were  obtained  in  June-July  1983  during  the 
MIZEX  Greenland  Sea  Summer  Pilot  Study  in  which  both  authors 
participated  on  Norweigian  research  vessel  Polarbjorn.  This 
experiment  is  refered  to  as  "MIZEX-83." 


2.  MARGINAL  ICE  ZONE  (MIZ) 


2.1.  General 

The  vicinity  of  the  open  ocean  gives  the  marginal  ice 
zone  its  peculiar  properties.  Heat  is  transported  to  the  MIZ 
from  the  warm  sea  through  currents,  the  drift  of  the  ice  itself 
aid  via  the  atmosphere.  Penetration  of  swell  into  the  MIZ  tends 
to  break  ice  floes  into  small  pieces.  Due  to  the  presence 
of  the  open  boundary  the  ice  has  more  freedom  of  motion  than 
in  the  central  pack.  Specifically  it  is  easy  to  diffuse  the 
ice  edge  since  pack  ice  has  a  very  small,  if  any,  tensile  strength. 

On  the  other  hand  both  the  atmosphere  and  the  ocean  feel  a 
sharp  change  in  their  boundary  surface  at  the  ice  edge  from  which 
arise  a  number  of  phenomena  typical  to  this  region.  The  characteristic 
scale  of  the  MIZ  width  is  100  km  which  is  also  representative 
to  atmospheric  and  oceanic  processes  at  the  ice  edge. 

a.  MIZ  ice  features 

Three  sub-zones  can  be  distinguished  in  the  MIZ  (Bauer 
and  Martin,  1980):  edge  zone,  transition  zone  and  interior 
zone  (Fig.  2).  This  structure  is  caused  by  the  gravity  waves 
in  the  ocean.  In  the  edge  zone  the  ice  meets  the  full  wave 
spectrum  and  ice  blocks  may  even  become  piled  up  on  each  other. 

Floes  are  rectangular  and  of  tens  of  meters  in  diameter  in 

the  transition  zone  where  waves  cause  floe  breakup  and  some 

rafting.  Ocean  swell  can  be  observed  in  the  interior  zone 

but  it  has  no  more  the  energy  to  fracture  ice  floes.  Thus 

waves  set  up  a  peculiar  floe  size  distribution  in  the  MIZ  although 

pack  ice  always  does  consist  of  a  number  of  finite  separate 

floes.  The  division  of  Bauer  and  Martin  was  based  on  the  Bering 

Sea  ice  conditions  but  seems  to  be  generally  valid  in  the  qualitative 

sense . 

The  mechanism  of  wave  induced  ice  breakup  in  the  MIZ  is 
explained  in  Wadhams  (1980  and  1984)  and  Squire  (1983)  .  As  a  normal  gravity 
wave  is  coupled  to  an  ice  sheet  it  propagates  partly  as  a  sinusoidal 
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flexural-gravity  wave  in  the  ice.  Fracturing  will  then  occur 
where  the  bending  stress  becomes  large  enough.  Short  wavelengths 
attenuate  much  more  rapidly  than  the  long  ones  and  only  the 
longest  ocean  swell  propagates  beyond  the  edge  zone.  Finite 
floes  have  a  rigid-body  response  to  a  wave  field  which  relieves 
the  bending  stress  on  a  floe,  and  for  a  given  wave  frequency 
the  wave  height  required  to  cause  fracturing  increases  as  the 
floe  diameter  decreases. 

The  wave- induced  breakup  would  thus  tend  to  produce  an 
increasing  floe  size  with  distance  from  the  edge  but  in  reality 
a  discontinuous  zone  structure  is  observed.  Wadhams  (1980) 
suggested  that  the  disagreement  is  due  to  mixing  of  floes  by 
fluctuating  ice  motion.  He  further  describes  ice  margin  as 
a  scrapyard  which  pulverizes  any  large  floes  which  venture 
into  its  sphere  of  influence.  Another  significant  effect  of 
the  ocean  waves  is  that  they  act  to  maintain  a  compact  ice 
edge  because  diffusion  of  the  edge  is  opposed  by  the  radiation 
pressure  of  the  wave  field  (Wadhams,  1980). 

In  the  MIZ  the  melting  of  ice  is  enhanced.  First,  there 
is  often  a  polar  front  in  the  ocean  near  the  ice  edge  and  the 
temperature  of  the  surface  water  may  change  by  several  degrees 
in  a  short  distance.  Eddies  develop  along  the  front  and  transport 
warm  water  into  the  sea  beneath  the  MIZ,  and  ice  may  be  driven 
across  the  front  by  the  eddies  or  an  off-ice  wind.  Secondly, 
the  presence  of  a  free  boundary  in  the  MIZ  allows  rapid  divergence 
under  off-ice  forcing.  Open  patches  and  leads  absorb  much 
more  solar  radiation  than  the  ice  and  the  heat  gain  is  used 
to  lateral  melting  of  ice  floes.  This  is  clearly  effective 
in  summer.  On  the  other  hand,  in  winter  openings  freeze  rapidly 
and  thus  increase  the  production  of  ice.  Wadhams  et  al.  (1979) 
stated  that  near  the  ice  margin  the  penetrating  swell  generates 
an  oscillating  shear  current  at  the  ice-water  interface  which 
further  increases  the  melt  rate  (when  the  water  temperature 
is  above  the  freezing  point) . 

The  decrease  of  ice  thickness  across  the  MIZ  was  clearly 


shown  by  Wadhams  et  al.  (1979)  (Fig.  3).  The  mean  ice  thickness 
was  doubled  at  the  distane  of  about  100  km  from  the  ice  edge. 

In  this  particular  location  the  thickness  differences  are  partly 
due  to  the  difference  in  ice  origin  (Wadhams,  1983b).  In  the 
Bering  Sea  MIZ  the  melt  rate  may  be  as  high  as  1  m  per  day 
as  the  ice  drifts  across  the  front  (Bauer  and  Martin,  1980; 

Josberger,  1983). 

Due  to  the  breakup  and  melting  of  ice  the  erosion  of  deformation 
features  such  as  pressure  ridges  is  pronounced  in  the  MIZ. 

The  observations  of  Kozo  and  Diachok  (1973)  in  the  Denmark 
Strait  show  a  rapid  decrease  in  the  frequency  and  size  of  ridges 
across  the  MIZ  (Figs.  4  and  5).  They  also  concluded  that  the 
sail  height  and  keel  depth  distributions  are  similar  and  scale 
by  a  constant  ratio  of  1:5  (Fig.  5).  As  stated  by  Wadhams 
(1980)  an  important  consequence  of  the  smoothing  of  roughness 
elements  is  that  the  form  drag  coefficients  of  wind  and  water 
become  reduced.  However,  the  net  effect  to  ice  drift  speed 
is  not  clear  since  the  free  drift  velocity  is  roughly  proportional 
to  the  ratio  of  the  wind  drag  coefficient  to  the  water  drag 
coefficient . 

Mixing  of  ice  floes  is  efficient  in  the  MIZ  due  to  the 
relatively  high  freedom  of  motion.  As  the  ice  cover  opens 
up,  ice  floes  tend  to  drift  as  individuals,  and  variations 
in  floe  form,  floe  mass  and  forcing  cause  large  differential 
velocities.  Giant  floes  and  icebergs  may  enter  the  MIZ  and 
have  a  yet  unknown  effect  on  the  drift  pattern. 

b.  Ice  motion  in  the  MIZ 

Towards  the  ice  edge  the  freedom  of  ice  movement  increases, 
and  higher  magnitude  and  variability  in  ice  velocity  is  likely 
to  occur.  Due  to  the  pressence  of  the  free  boundary  ice  compactness 
may  vary  much  which  means  pronounced  variations  in  ice  strength 
since  the  logarithm  of  that  is  proportional  to  the  compactness 
(Hibler,  1979).  In  the  case  of  very  low  strength  the  ice  drift 
follows  the  local  external  forcing.  With  increasing  strength 
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the  drift  follows  the  forcing  integral  over  a  scale  increasing 
up  to  1()3  km  in  the  Central  Arctic. 

First  a  few  observed  cases  are  considered.  A  good  illustration 
of  the  velocity  increase  across  the  MIZ  has  been  presented 
by  Vinje  (1977)  on  the  basis  of  Landsat  images  (Fig.  6).  He 
gave  a  wind-driven  and  an  ocean-driven  displacement  over  about 
two  weeks  and  a  general  increase  over  the  MIZ  is  seen  in  both 
cases.  Johannessen  et  al .  (1983b)  showed  four  cases  of  daily 
average  MIZ  flow  where  there  was  a  very  sharp  increase  in  velocity 
within  about  5  km  from  the  ice  edge  (Fig.  7).  Its  cause  was 
uncertain.  Lepparanta  (1982)  gave  a  daily  deformation  of  the 
Baltic  Sea  ice  cover  and  the  magnitude  increased  by  a  factor 
of  2  to  3  near  the  ice  edge  (Fig.  8).  In  this  case  it  was 
seen  that  the  ice  pack  got  more  open  when  approaching  the  edge. 

There  are  certain  features  in  the  MIZ  which  seem  to  have 
a  significant  influence  on  its  dynamics.  They  are:  1)  Small 
floe  size;  2)  Decrease  of  floe  size  when  approaching  the  edge; 

3)  Decrease  of  ice  thickness  when  approaching  the  edge;  4) 

Presence  of  free  boundary;  5)  Coupling  with  atmosphere  and 
ocean  which  feel  a  sharp  change  in  their  boundary  surface; 

6)  Penetration  of  ocean  swell;  7)  Smoothing  of  roughness  elements 
when  approaching  the  edge.  It  is  a  difficult  task  for  the 
MIZEX-program  to  try  to  solve  the  role  of  these  factors. 

Presently  a  number  of  qualitative  characteristics  of  the  MIZ 
flow  are  known  and  their  causes  are  widely  discussed.  Hence 
it  is  considered  important  to  list  and  explain  them  here: 

1)  Edge  dispersion.  It  has  been  noted  often  that  once  the 
forcing  turns  off-ice,  the  ice  edge  diffuses. 

2)  Banding.  It  may  happen  that  narrow  bands  form  perpendicular 
to  the  wind  direction  when  the  wind  is  off-ice.  These  bands 
seem  to  conserve  their  form  remarkably  well.  (Zubov,  1945; 

Martin  et  al.,  1983). 

3)  Edge  compaction.  For  on-ice  forcing  ice  floes  drift  together 
and  the  compactness  practically  drops  from  one  to  zero  at  the 
ice  edge.  This  feature  is  quite  frequently  noted. 


A)  Little  ridging.  It  seems  that  new  ridges  are  not  or  very 
seldom  formed  in  the  MIZ  except  for  some  piling-up  of  ice  blocks 
due  to  ocean  waves  in  the  edge  zone. 

5)  Ice  edge  jet.  There  is  some  indication  that  the  ice  velocity 
may  sharply  increase  near  the  ice  edge.  However,  it  is  not 
certain  what  is  the  main  mechanism  which  may  cause  the  jet. 
(Johannessen  et  al.,  1983b;  Roed  and  O'Brien,  1981;  Lepparanta, 
1984;  Lepparanta  and  Hibler,  1984b). 

6)  Forcing  of  upwelling  at  the  ice  edge.  When  the  surface 
stress  on  the  ocean  is  parallel  to  the  ice  edge,  a  discontinuity 
in  it  may  occur  at  the  edge  and  cause  upwelling.  It  is  an 

ice  dynamics  problem  to  solve  whether  the  discontinuity  exists 
and  what  is  its  sign.  (Buckley  et  al.,  1979;  Clarke,  1978; 

Roed  and  O'Brien,  1983). 

The  first  and  second  points  are  concerned  with  a  diverging 
ice  edge.  The  ice  offers  then  little  resistance  and  is  more 
or  less  a  passive  tracer  of  the  forcing.  Thus  these  cases 
are  likely  to  be  externally  caused.  The  third  and  fourth  points, 
on  the  other  hand,  are  clearly  a  manifestation  of  ice  interaction. 
The  stress  within  the  MIZ  does  not  get  very  high  because  of 
the  free  boundary,  and  the  compressive  strength  of  ice  prohibits 
further  convergence  once  the  compactness  has  come  close  enough 
to  one.  The  jet  is  in  principle  possible  due  to  ice  interaction, 
external  forcing  or  both.  Upwelling  at  the  ice  edge  needs 
a  discontinuity  and  hence  a  compact  ice  edge,  which  means  that 
the  ice  interaction  then  has  a  significant  role  in  ice  motion. 

It  seems  likely  that  the  wind  stress  is  greater  over  the  ice 
than  over  open  water.  Due  to  frictional  losses  within  the 
ice  only  a  part  of  the  wind  stress  is  transmitted  through  the 
ice  cover  into  the  ocean.  Thus  either  sign  of  the  ocean  forcing 
discontinuity  is  possible  at  the  ice  edge. 

c.  MIZ  ice  rheology 

Study  of  ice  rheology  is  involved  with  the  nature  of  ice 
interaction  and  aims  at  establishing  a  law  for  determining  the 
stress  within  the  ice  pack.  It  can  be  generally  stated  that  the 
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ice  stress  depends  on  the  strain  rate  and  thickness  distribution 
(see,  e.g.  Hibler,  1979).  The  dynamical  meaning  of  the  rhe¬ 
ological  law  is  that  it  tells  how  kinetic  energy  is  dissipated 
and  transmitted  within  the  ice.  Formally,  the  divergence  of 
the  ice  stress  appears  as  a  term  in  the  momentun  equation. 

There  are  three  classical  pack  ice  rheologies:  1)  Free 
drift  (Nansen,  1902);  2)  Linear  viscosity  (Laikhtman,  1958;  Glen, 

1970);  3)  Plasticity  (Coon  et  al . ,  1974;  Hibler,  1979).  A  variation 
of  these  cases  is  a  nonlinear  viscous  model,  which  is  essentially 
the  approach  of  Hibler  (1979).  They  are  illustrated  in  Fig.  9. 

The  free  drift  assumption  simply  states  that  there  are  no 
stresses  within  the  ice  and  is  valid  once  the  ice  pack  opens  up 
so  that  floes  do  not  interact  with  each  other.  Such  is  the  case 
sometimes  in  summer  conditions  as  has  been  recently  shown  by  McPhee 
(1980).  Once  the  ice  pack  becomes  compact  the  free  drift  approximation 
fails  since  it  does  not  give  any  resistance  to  deformation.  E.g., 
an  infinite  build-up  of  ice  would  occur  at  the  coastline  during 
an  on-shore  wind. 

Linear  shear  viscous  rheology  was  used  by  Laikhtman  (1958)  to 
describe  the  flow  of  ice  near  a  fixed  boundary  and  Campbell  (1965) 
could  simulate  with  it  the  general  features  of  the  ice  drift  in  the 
Arctic  Ocean.  This  original  linear  model  considered  only  the  shear 
viscosity  and  Glen  (1970)  extended  it  to  include  the  bulk  viscosity 
as  well.  The  linear  viscous  approach  is  still  being  applied  to 
some  degree  due  to  its  computational  easiness.  Lepparanta  (1981a) 
has  used  Glen's  extended  law  in  a  routine  ice  forecasting  model 
for  the  Baltic  Sea  with  some  success. 

In  the  1970s  the  linear  viscous  law  was  highly  criticized. 

The  existence  of  a  significant  yield  strength  seemed  to  be  clear 
and  the  coastal  boundary  layer  of  the  ice  pack  highly  nonlinear. 

The  latter  feature  appears  as  the  well-known  shear  zone,  e.g., 
off  the  northern  coast  of  Alaska  and  Greenland  in  the  Arctic. 
Additionally,  ridging  is  considered  as  the  dominant  dissipation 
mechanism  in  the  ice  (Rothrock,  1975)  and  in  the  small  scale  ridging 
has  a  plastic  nature  (Parmerter  and  Coon,  1972).  Consequently 
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a  plastic  rheology  was  suggested  (Coon  et  al.,  1974).  Due  to 
mainly  computational  reasons  the  plastic  ice  rheology  has  been 
approximated  by  an  elastic-plastic  law  (Coon  et  al.,  1974)  and 
a  viscous  plastic  law  (Hibler,  1979).  Both  have  led  to  very 
good  results  in  the  Arctic  Basin  (Hibler,  1979;  Pritchard,  1980). 

An  additional  important  feature  of  these  plastic  (or  nonlinear 
viscous)  rheologies  is  that  they  include  a  thickness-dependent 
ice  strength  and  thus  describe  pack  ice  as  a  work-hardening 
material.  This  gives  a  strong  restriction  to  the  amount  of  ice 
build-up.  Hibler  (1980b)  has  got  a  very  good  agreement  between 
his  model  thicknesses  and  observations. 

In  the  MIZ  ice  properties  and  processes  have  specific  chara¬ 
cteristic  features,  and  it  is  not  clear  which  would  be  a  proper 
rheological  model  for  it.  The  relevant  question  is  that  what 
is  the  main  dissipation  mechanism  in  the  MIZ.  Since  ridging 
seems  to  be  ruled  out,  likely  candidates  should  be  friction 
between  floes  in  contact  and  floe  collisions.  This  has  very 
interesting  implications.  In  the  general  theory  of  flow  of 
granular  media  the  former  case  is  sometimes  named  slow  and  the 
latter  rapid  flow.  It  has  been  shown  that  the  slow  granular 
flow  is  plastic  (e.g. ,  Mandl  and  Luque ,  1970)  whereas  in  the 
rapid  flow  the  stress  depends  quadratically  on  the  strain  rate 
(Shen  and  Ackermann,  1982).  The  difference  between  these  two 
laws  is  illustrated  in  Fig.  10.  It  is  possible  that  the  MIZ 
flow  may  obey  either  of  them  depending  on  the  flow  characteristics, 
or  even  something  of  mixed  type. 

2.2  Summer  MIZ  in  the  Fram  Strait 
a.  General  description 

On  average  at  the  end  of  June  the  coast  of  Svalbard  is 
open  west  of  20°E  (Fig.  11).  The  annual  variation  of  the  mean 
monthly  position  of  the  ice  edge  is  about  200  km.  Short  term 
variation  can  be  comparable  to  that  since  a  lOm/s-wind  may  displace 
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ice  by  20  km  in  a  day. 

An  old  description  of  Koch  (1945)  says  that  the  ice  in 
the  eastern  part  of  the  Fram  Strait  originates  from  the  Eurasian 
side  of  the  Arctic  Ocean  and  is  not  very  heavy.  The  western 
part,  on  the  other  hand,  contains  old  and  thick  ice  from  the 
Central  and  Western  Arctic.  The  idea  is  based  upon  assuming 
a  laminar  ice  outflow  from  the  Arctic  Ocean.  Satellite  images 
have  shown  large  fluctuations  in  ice  motion  which  indicates 
efficient  mixing  of  ice  floes  but,  however,  a  recent  work  of 
Wadhams  (1983b)  supports  Koch’s  description  to  some  degree. 

Wadhams  (1983b)  presented  ice  thickness  distribution  for 
the  Fram  Strait  showing  gradients  steeper  than  1  m  per  100  km 
(Fig.  12).  Our  observations  at  MIZEX-83  drift  phase  location 
fit  well  in  his  map.  In  the  summer  of  1981  Vinje  et  al.  (1983) 
found  heavily  ridged  ice  close  to  the  ice  edge  but  during  the 
proceeding  months  there  had  been  an  unusual  east-southeasterly 
drift  of  ice  north  of  Fram  Strait. 

Fram  Strait  is  the  main  channel  of  the  outflow  of  ice  from 
the  Arctic  Ocean.  It  has  been  estimated  that  the  average  value 
is  about  10^  m3  s“l  (Vinje,  1976).  As  the  ice  leaves  the  Arctic 
Ocean  it  enters  the  East  Greenland  current  and  its  speed  increases 
from  6-12  cm/s  at  80-81°N  to  6-37  cm/s  at  76-79°N  (Vinje,  1976). 
The  ice  velocity  tends  to  increase  also  across  the  Fram  Strait 
so  that  the  maximum  value  occurs  at  the  ice  edge  (see  Fig.  6). 

For  a  more  complete  discussion  of  the  ice  in  the  northern 
Greenland  Sea  the  reader  is  referred  to  Wadhams  (1980),  Weeks 
(1982)  and  Overgard  et  al.  (1983). 

b.  Conditions  during  the  MIZEX-83  experiment 

In  the  summer  of  1983  the  ice  conditions  seemed  to  be  normal 
in  the  northern  Greenland  Sea  (Fig.  13).  The  present  authors 
participated  in  the  MIZEX-83  experiment  on  R/V  Polarbjorn  for 
carrying  through  ice  kinemetics  studies.  These  were  done 
during  the  so  called  drift  phase,  June 


28  to  July  9,  when  Polarbjorn  was  moored  to  an  ice  floe  in  the 
northern  part  of  the  Fram  Strait  (Fig.  13).  However,  before 
the  drift  phase  began  the  ship  sailed  in  the  MIZ  for  several 
days  in  order  to  deploy  instruments  on  and  beneath  the  ice. 

Thus  notes  on  the  features  of  the  MIZ  could  be  made. 

Due  to  difficulties  in  navigating  in  the  MIZ  Polarbjorn 
could  not  go  deeper  than  20-25  km  into  the  MIZ.  The  ice  edge 
was  usually  well  defined  and  the  edge  zone  well  developed. 

The  ice  was  not  uniformly  distributed  but  clustering  of  floes 
was  noted  as  a  typical  phenomenon.  In  places  there  were  open 
patches  and  in  places  closely  packed  clumps  of  floes.  This 
is  probably  important  to  the  initialization  of  edge  dispersion. 

As  an  example,  visual  observations  during  one  transect  are  illustrated 
in  Fig.  14.  During  the  drift  phase  Polarbjorn  was  located  at 
about  20  km  from  the  edge  and  the  ice  was  there  more  uniformly 
distributed  with  compactness  greater  than  about  0.7. 

Ice  floes  were  typically  10  to  100  meters  in  diameter  and 
there  was  some  average  size  increase  with  the  distance  from 
the  edge.  Deeper  in  the  pack  large  floes  appeared.  In  the  drift 
phase  study  area  the  floe  size  was  quite  homogeneous,  around 
50  meters. 

The  MIZ  ice  was  a  mixture  of  first-year  and  multi-year 
ice  and  the  ice  thickness  ranged  from  1  to  4  meters.  It  seemed 
that  the  amount  of  first-year  ice  was  greater.  Quite  high, 
even  3  meters,  ice  ridge  sails  were  seen.  One  could  definitely 
say  that  the  ice  was  much  more  heavily  ridged  than  reported 
by  Overgaard  et  al.  (1983)  for  north  of  Nordaustlandet  in  the 
summer  of  1980. 

On  one  occasion  the  efficiency  of  floe  mixing  was  noted. 

Two  marked  ice  floes,  originally  side  by  side,  were  found  2 
km  apart  from  each  other  after  three  days. 
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3.  Modeling  ice  dynamics  in  the  MZZ 
3.1.  System  of  equations 

The  dynamics  of  pack  ice  is  governed  by  the  equations  of  con¬ 
servation  of  momentum  and  mass  (e.g.  Hibler,  1984) : 

+  +  7„  +  'a  +  ?g  ' 

§|  =  -GV-u  +  , 

where  is  the  density  of  ice,  h  mean  ice  thickness,  D/Dt 
material  time  derivative,  u=(u,v)  ice  velocity,  f  Coriolis  para¬ 
meter,  k  the  unit  vector  vertically  upward,  a  ice  stress  tensor, 
Ta  w*nc*  stress»  tw  water  stress,  gravitational  force  due  to 
the  sea  surface  tilt,  G  ice  thickness  distribution  function  and 
<|>  mechanical  thickness  redistributor.  Various  forms  of  Eq.  (3.1) 
have  been  used  for  a  long  time.  The  continuity  equation  above 
was  presented  by  Thorndike  et  al.(1975);  G(h)  is  the  area  of  ice 
thinner  than  h  and  the  function  Q  describes  the  transfer  of  ice 
from  one  thickness  category  to  another  due  to  deformation. 

The  main  ice  mechanics  problems  are  to  establish  a  proper 
formulation  for  the  ice  stress  and  the  mechanical  redistributor. 
In  general  one  has 

a  =  a (G, e )  ,  (3.3) 

♦  “  'I'  (G,  e)  ,  (3.4) 

where  e  =  1/2 [Vu  +  (Vu)^]  is  the  strain  rate  tensor  of  the  ice. 
Once  explicit  expressions  are  given,  the  system  becomes  diff¬ 
icult  to  solve.  Eqs/3.1-2)  are  coupled  since  the  ice  stress 
depends  on  the  thickness  distribution.  Note  that  there  is  also 
a  weak  coupling  even  in  the  absence  of  ice  stress  through  the 
appearance  of  the  mean  ice  thickness  in  Eq. (3.1). 

The  forcing  terms  are  expressed  as 


(3.1) 

(3.2) 


(3.5) 


T  =  P  C  | U  |  (cos$  U-  +  sin$  k*U_}  , 

a  a  a  a  a  a 

7W  =  PWCW | Uw-u! (cose  (Uw-u)  +  Sine  kx(Uw-u)>  ,  (3.6) 

Fg  =  -pj^hg^C  ,  (3.7) 

where  t  and  tw  are  the  air  and  water  densities,  and  Cw  air 
and  water  drag  coefficients,  Ua  and  Uw  air  and  water  velocities, 

<(>  and  0  boundary  layer  turning  angles  for  air  and  water,  g  is  the 
gravity  acceleration  and  5  the  dynamic  height  of  the  sea  surface. 
In  analytic  studies  it  is  often  convenient  to  linearize  the 
water  stress  using  the  Ekman  boundary  layer  solution  modified 
by  the  angle  of  ir/4  -  0, 

Tw  =  pw^  (cos6  (Uw~u)  +  sin0  icx  (Uw~u)  )  ,  (3.8) 

where  <  is  the  kinematic  eddy  viscosity  coefficient  for  water. 

When  modeling  MIZ  ice  dynamics  two  meaningful  simplifications 
can  be  made.  First,  the  MIZ  is  idealized  as  a  long  narrow  zone 
with  straight  and  parallel  boundaries  (Fig. 15),  and  it  is 
assumed  that  spatial  variations  are  much  larger  perpendicular 
than  parallel  to  the  ice  edge.  Then,  aligning  the  y-axis  with 
the  ice  edge (Fig. 15)  we  have 


This  is  a  similar  boundary-layer  model  to  Laikhtman's  (1958)  and 
has  been  applied  for  the  MIZ  by  Roed  and  O'Brien  (1981)  and  Lep- 
paranta  and  Hibler  (1984b). 

Secondly,  the  thickness  distribution  can  be  reduced  to  two 
levels:  the  fraction  of  open  water  and  the  fraction  of  ice  cover 
represented  with  one  thickness  value  (in  each  element).  This  can 
be  reasoned  by  the  fact  that  we  consider  dynamics  only  and  need 
not  worry  about  the  different  growth  rates  of  thin  ice  and  thick 
ice.  In  addition,  the  changes  in  the  thickness  distribution  due 
to  deformation  are  caused  by  opening  and  closing  of  leads  more 


than  by  ridging  in  the  MIZ.  The  continuity  equation  in  the 
two-level  case  can  be  easily  transformed  into  one  for  the  mean 
thickness  h  and  one  for  the  ice  compactness  A  as  has  been  done 
by  Doronin  (1970).  Henceworth  the  tilde  is  omitted  from  the 
notation  of  mean  thickness.  The  role  of  the  redistributor  is 
now  to  take  care  that  the  conditions  h>0  and  0<A^1  are  always 
fullf illed. 

Finally,  we  shall  consider  time-scales  more  than  one  hour  or 
so.  The  inertial  term  is  known  to  be  small  for  such  conditions 
and  is  hence  neglected.  Then  the  momentum  equation  has  only  one 
independent  variable  left.  Our  MIZ  ice  dynamics  model  has  now 
become 


3X  ^°XX 

,  c  )  - 

'  *y 

Pihfkxu  +  tw  +  xa  +  Fg  =  0  , 

(3.10) 

3h  _ 

3 1 

h|a  ,  h>_o. 

(3.11) 

3A 

3t 

-St- 

Af^ 

(3.12) 

Additionally,  the  assumption  (3.9)  simplifies  the  strain  rate 
tensor  to  the  form 


e  = 


3u 

3x 


1 

2 


3  V 
dx 


(3.13) 


As. the  initial  conditions  the  distribution  of  ice  compactness 
and  thickness  across  the  MIZ  must  be  spesified.  The  ice  edge 
location,  n=Ut),  is  initially  at  x=0.  The  boundary  between  the 
MIZ  and  fast  ice  (or  interior  ice)  is  fixed  at  x=L  and  the 
velocity  of  ice  is  there  given,  u=uL,  i.e.  in  the  case  of  fast 
ice  u^=0*  At  the  ice  edge  the  ice  stress,  thickness  and  compact¬ 
ness  vanish.  Formally, 

u  =  uL  , 

o=0,  h  =  0,  A  =  0. 


x=L: 
X*  lS 
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(3. 14. a) 
( 3. 14. b) 


The  use  of  a  fixed  boundary  is  reasonable  when  there  is  land  or 
fast  ice  within  a  few  hundred  kilometers  from  the  ice  edge  as  in 
the  Greenland  Sea.  On  the  other  hand  if  such  is  not  the  case  the 
model  represents  a  two- layer  approximation  to  the  MIZ  -  interior 
ice  system  where  stress  is  allowed  to  occur  at  the  layer  bound¬ 
ary.  The  two-layer  structure  might  be  caused  by  a  rapid  change 
in  ice  characteristics  or  forcing.  An  other  possibility  would  be 
to  replace  the  condition  (3. 14. a)  by  an  infinite  boundary  con¬ 
dition. 

The  no  stress  condition  at  the  ice  edge  might  be  replaced  by 
equating  the  normal  stress  to  the  radiation  pressure  of  the  ocean 
waves.  The  directional  spectrum  of  the  wave  field  just  outside 
the  ice  edge  is  directed  towards  the  ice,  and  the  pressure  force 
per  unit  length  on  the  ice  edge  is  (Longuet-Higgins,  1977; 
Wadhams,  1980) 


Fw  "  T  pw*<a2  "  ar  +  at> 


(3.15) 


where  a,  ar  and  afc  are  the  incident,  reflected  and  transmitted 


wave  amplitudes,  respectively.  This  force  varies  between  zero 


and 


pwga  /2’ 


Let  us  now  write  Eq. (3.10)  in  component  form: 


ax  °xx  +  Pihfv  “  Pwcwluw_ul  *CO80  u  ~  sine  v) 


(3. 16. a) 


+  pwCwlUw-ul  <cos9  Uw”  sine  V  +  Tax  +  Fgx  =  0  ' 


ax  axy  "  Pihfu  ~  pwCwlUw'ui  (cose  v  +  sine  u> 


+  pwCw|Uw_u| (cos9  Vw  *  sin9  Uw}  +  Tay  =  0 


(3. 16. b) 


Note  that  the  linear  form  (Eq.3.8)  for  the  water  stress  would 
separate  it  into  independent  damping  and  forcing  terms.  The 
solution  of  Eq. (3.16)  is  the  steady  state  solution  of  ice 
velocity  for  a  given  mass  distribution  and  is  valid  in  time 
scales  from  one  hour  up  to  perhaps  a  day  depending  on  the  con- 


n .  . .  ^ 

ditions  in  each  case.  In  longer  time  scales  the  adjustment 
between  the  mass  and  velocity  distribution  may  significantly 
modify  the  flow.  To  obtain  the  steady  state  solution  for  the 
adjusted  case  Eq. (3.16)  must  be  satisfied  together  with 


3h  „ 

at  =  0 


3A 

3t 


=  0 


(3.17) 


It  follows  directly  from  Eqs. (3. 11-12)  and  from  the  presence  of 
the  ice  edge  that  the  condition  (3.17)  is  equivalent  to 


u  =  0  . 


(3.18) 


Thus  to  solve  the  adjusted  case  we  only  need  to  delete  u-velocity 
in  Eq. (3.16) . 

The  standard  parameters  for  the  modeling  studies  are  shown  in 
Table  1.  The  atmospheric  and  oceanic  parameters  correspond  to 
geostrophic  motion  beyond  the  layers  of  frictional  influence  of 
the  ice.  The  eddy  viscosity  coefficient  for  water  has  been 
selected  so  that  the  linear  and  quadratic  stress  laws  coincide 
at  |uw~u|  =10  cm/s. 


Table  1.  Standard  modeling  parameters. 


g  =  9.8  ms 


-2 


f  =  1. 46*10 


-4  -1 


Cw  «  5.5x10 
0  »  25° 


-3 


C  >  1.2x10 

cl 

*  =  25° 


-3 


-3 


=  910  kg  m 
pw  =  103  kg  m-3 


3.2.  Free  drift 


=  21  cm  s 


-2 


p  _  =  1.3  kg  m 

cl 


-3 


The  steady  state  free  drift  solution  of  the  momentum  equation 
is  quite  straightforward.  The  equation  reads  as 


-pjhfkxu  +  t  +  t  +  f, 
i  awe 


=  0  . 


(3.19) 
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Employing  the  geos trophic  balance  m  the  ocean  the  tilt  term  can 
be  expressed  as  F^  =  P^hfkxUw.  Now,  it  is  easy  to  see  (e.g. 
Hibler,  1984)  that  the  free  drift  velocity  is  obtained  from 


-Pihfkxur  +  Tw(ur)  +  ta  =  0  , 


u  =  ur  +  Uw  ' 


(3. 20. a) 
(3. 20.b) 


where  ur  is  the  velocity  of  ice  relative  to  the  geostrophic 


current  and  the  notation  Tw(ur)  means  that  u-Uw  is  replaced  by 


ur  in  the  water  stress  expression. 

First  we  take  the  linear  water  stress  law  (Eq.3.8).  Eq. (3. 20. a) 
becomes  in  matrix  form 


/  cose 

-sine\  , 

f  0 

\  ^ 

-1 

-  pw^ 

w  l  sine 

^  ' 

W  p±hf  I 

cose  / 

1 

0  J 

(3.21) 


This  can  be  written  as 

(  C1  0  ^  ( cosei 

-sin^  \ 

\  0  c1y^sin81 

cosB^  J 

tos  * 

p.C  |U  1 

|  1 

a  a 1  a 

Is  in* 

where 

6^  =  arccos{ [1 

+  (tanO 

/  u  x 

!  r 


r  ; 


(3.22) 


-sm 
cos*  j  y 


v_  ! 


Pihf 


■>2r1/2>  , 


(3. 23. a) 


COS0 

w 


c,  =  p../7f  cose  cosB,  . 
1  w  1 


(3. 23. b) 


Through  inverting  the  matrixes  on  the  left  hand  side  of  Eq. 

(3.22)  we  see  that  in  the  linear  free  drift  solution  the  wind 
velocity  is  contracted  by  paCalua|/c^  and  rotated  by  *  ~  6^. 

The  case  of  quadratic  water  stress  is  more  difficult  and 
does  not  give  an  explicit  solution.  Lepparanta  (1981b)  multiplied 


Eq. (3. 20. a)  scalarly  with  itself  and  with  ur,  which  led  to  a 
pair  of  nonlinear  equations  for  the  dimensionless  speed  |ur|/|Ua| 
and  the  deviations  angle  between  the  ice  and  wind  velocities,  6. 
The  solution  depends  on  the  boundary  layer  angles  and  two 
dimensionless  products  which  were  chosen  a? 


»  -  i 

Hw  w 


R  = 


Pjhf 


iu  | /p  c  p  c 

i  a1  Ma  aMw  w 


Then, 


lurl/lual  ■  Na*  F (R, 0 ) 


cos(<)>+8)  =  cosQ-F  (R,0)  , 


where  the  function  F  is  defined  by 


F2  +  2sin8  RF3  +  R2F2  -  1 


=  0 


(3.24) 


(3. 25. a) 
(3. 25. b) 


(3.26) 


For  a  fixed  oceanic  boundary  layer  angle  0,  F  decreases  monoton- 
ically  from  1  to  0  as  R  goes  from  zero  to  infinity  (Fig. 16).  Note 
that  F  is  only  weakly  dependent  on  9. 

The  basic  qualitative  difference  between  the  linear  and  qua¬ 
dratic  cases  is  in  the  wind  dependency  of  the  solution.  The  linear 
case  states  that  the  ice  speed  is  proportional  to  the  squared 
wind  speed  whereas  in  the  quadratic  case  ice  speed  increases  only 
slightly  stronger  than  in  direct  proportion  to  the  wind  speed. 

The  deviation  angle  depends  naturally  on  the  atmospheric  boundary 
layer  angle  but  in  the  quadratic  case  it  also  decreases  with  in¬ 
creasing  wind  stress.  Both  cases  give  qualitatively  similar  role 
to  the  ice  thickness  h:  as  h  increases,  the  ratio  of  ice  speed 
to  wind  speed  decreases  and  the  deviation  angle  increases. 

The  free  drift  approximation  is  reasonable  when  the  compactness 
of  ice  is  less  than  about  0.8  and  it  is  hence  useful  for  studying 
the  edge  dispersion.  Even  this  simple  law  leaves  much  freedom  to 
variability  of  ice  motion  through  forcing,  boundary  layer  param¬ 
eters  and  the  thickness  of  ice.  Here  we  concentrate  on  the  last 
of  these  factors  since  its  variability  is  rather  well  known;  e.g. 
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during  MIZEX-83  it  was  noted  that  the  Greenland  Sea  MIZ  is  a 
mixture  of  first-year  and  multiyear  ice  floes  with  thicknesses 
ranging  from  1  to  4  m.  The  basic  idea  of  the  thickness  effect 
dates  back  to  Zubov (1945)  and  has  been  recently  applied  for  the 
MIZ  by  Lepparanta  and  Hibler  (1984a). 

The  thickness  variations  may  have  three  dynamically  important 
consequences:  1)  Mixing  of  ice  floes  as  they  drift  with  diff¬ 
erent  speeds;  2)  Clustering  of  floes  as  light  floes  catch  up 
with  heavy  floes  in  some  places  and  drift  apart  from  them  in 
others?  3)  Dispersion  as  light  floes  drift  fast  away  from  the 
MIZ.  It  is  anticipated  that  a  variable  form  drag  of  floes  might 
lead  to  qualitatively  similar  results  but  it  is  not  known  how 
the  form  drag  varies. 

To  understand  the  thickness  effect  quantitatively  some  free 
drift  solutions  are  shown  in  Fig. 17  in  vector  form.  Since  1  cm/s 
approximately  equals  1  km/day  it  is  seen  that  neighboring  first- 
year  and  multiyear  ice  floes  can  drift  apart  several  kilometers 
in  one  day  due  to  the  thickness  difference  only.  Such  a 
phenomenon  was  in  fact  observed  in  MIZEX-83.  In  particular, 
after  several  days  of  off-ice  winds,  two  originally  adjacent 
floes  that  had  been  marked  were  relocated  and  found  to  be  about 
2  km  apart. 

It  is  also  notable  that  floe  clustering  due  to  differential 
drift  of  floes  is  a  substantial  conceptual  change  from  the 
thickness  distribution  concepts  used  in  most  sea-ice  models 
(see  e.g.  Hibler,  1984).  In  conventional  sea-ice  models,  one 
velocity  is  normally  used  for  a  variety  of  ice  thicknesses  in  a 
region,  and  statistical  ice  thickness  variations  are  then  deduced 
from  the  gradients  of  this  velocity  field.  Taking  into  account 
differential  drift  of  different  floes  would  require  a  conceptual 
reformulation  of  such  theories. 

To  examine  whether  in  fact  such  drift  variations  could  create 
clusters,  a  Monte  Carlo  simulation  was  carried  out.  The  results 
clearly  support  our  reasoning  above  (Fig. 18).  Initially,  ice 
floes  were  randomly  distributed  in  space  with  random  thicknesses 
up  to  5  m.  Only  one  drift  component  was  simulated.  The  system 
works  so  that  when  a  light  floe  catches  a  heavy  floe,  they  stick 


together  and  drift  as  a  group;  their  mean  thickness  then 
determines  the  drift  velocity.  Floe  clustering  is  seen  for  all 
our  ice  compactness  values.  Very  similar  results  were  obtained 
when  the  floes  were  initially  equally  rather  than  randomly 
spaced. 

There  is  not  yet  sufficient  data  on  the  drift  of  various  floes 
to  test  our  theory,  it  is  consequently  recommended  that  field 
measurement  relevant  to  this  problem  be  made  during  the  MIZEX-84 
main  experiment. 

Zubov  (1945)  suggested  that  the  thickness  variations  cause 
ice  bands.  This  may  not  be  exactly  true,  but  possibly  in  the 
initialization  of  band  formation  they  are  important.  Floe 
clustering  creates  openings  and  Wadhams  (1983b)  has  suggested 
that  the  wave  radiation  pressure  during  off-ice  winds  may  form 
bands  if  there  are  openings  present  near  the  ice  edge. 


3.3.  Analytic  linear  viscous  model 

The  steady  state  momentum  equation  of  pack  ice  with  linear 
viscous  rheology  is  now  applied  to  the  MIZ  model.  This  approach 
gives  a  possibility  to  obtain  analytical  solutions  for  certain 
interesting  cases.  There  has  been  much  speculation  on  the  MIZ 
ice  flow  on  the  basis  of  the  free  drift  law,  and  linear  models 
give  it  a  first-order  correction  which  considers  ice  interaction. 
In  the  present  cases  there  is  no  coupling  between  ice  thickness 
and  strength.  Consequently,  the  adjusted  steady  state  solution 
makes  sense  only  when  the  forcing  is  forced  so  that  u-velocity 
becomes  identically  zero.  This  is  a  far  too  spesific  restriction 
and  hence  the  adjusted  case  is  not  considered.  The  results  of 
this  section  are  based  on  a  recent  paper  of  Lepparanta  (1984). 

The  linear  viscous  rheology  is  expressed  as  (Glen,  1970) 

a  =  2ne  +  U~n)tre  I  ,  (3.27) 

where  I  is  the  unit  tensor,  and  n  and  C  are  the  shear  and  bulk 
viscosities  of  ice.  With  zero  y-derivatives , 


21 


d  ,  .  , du  dv, 

va  -  ^[(?+n)^  ,  n^]- 


(3. 28) 


We  take  the  linear  water  stress  formula  (Eq.3.8)  and  denote 


cw  =  ' 


(3.29 .a) 


F  =  t  +  c„(cos6  U  +  sine  k*U  )  +  F  ,  (3. 29. a) 

a  w  w  w  g 


Utilizing  Eqs. (3. 28-29) ,  the  momentum  equation  (Eq.3.10)  can  be 
written  in  component  form  as 


d  r„dvi 

d£[n35] 


c,  cose  u  + 
w 


c„cos8  v  - 
w 


(c  tsin9  + 
w 

(c  sin6  + 
w 


c  ^hf ) v  +  Fx 
Pihf ) u  +  Fy 


0  ,  (3. 30. a) 

0  .  ( 3. 30 . b) 


The  boundary  conditions  are  now  (from  3.14) 


x  =  0: 


..  \du 

(S+rr)^ 


(3. 31. a) 


x  =  L: 


u 


UL  '  v  =  VL  ’ 


(3. 31. b) 


a.  General  solution  with  constant  viscosities 

First,  to  shorten  our  notations,  Eqs. (3. 30.a-b)  are  divided 
by  the  constants  c+n  and  n  ,  respectively,  and  rewritten  as 


d2u 

— =-  +  a.  u  +  v.v  +  y . 
dx^  1  -1  L 


=  0  , 


(3. 32. a) 


d2v 

— 2  -  a2v  "  v2u  +  y2 
dx 


=  0  , 


(3.  32. b) 


where 
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t 


(c+nJa^  =  na2  =  CWCOS0 
(C+n)^  =  nv2  =  cwsin6  +  P^hf  ,  (3.33) 

U+n^  =  Fx,  nv2  =  Fy  . 

Using  the  elimination  method,  Eq.  (3.32)  can  be  combined  to  one 
for  u: 


4  2 

— j  -  (o^  +  +  ^ala2  +  viv2^u  ”  ^v1Y2  +  v2y1*  =  0 


(3.34) 


Once  this  is  solved,  v  is  obtained  directly  from  Eq. (3. 32. a). 
The  characteristic  polynomial  of  Eq. (3.34)  has  the  roots 


X,  0  ^  .  =  +  { 
1/  2/ 3, 4  — 


al+a2  .  r/al_a2v2  a/2,1/2 

— 2 —  +  [( — 2 — )  “  V1V2'  '  (3.35) 


and  thus  the  solution  is 


4 

u  =  £  C  exp(X.x)  +  u.  ,  (3.36) 

1=1  1  1 

where  u^  is  a  particular  integral  of  Eq.  (3.34)  and  Ce's  are 
undetermined  constants.  Then, 


v  = 


-1 


4 

£ 

e=l 


(oi_ 


X^C^expU^x) 


alul 


d2u 

dx2 


-  Y- 


}.  (3.37) 


Clearly,  the  roots  X ^ consist  of  4  real  values  when  the  expression 
in  brackets  in  Eq. (3.35)  is  non-negative;  this  is  true  when 


1  e 

2  /1+1/e 


p  ±hf 

_>  tan0  +  - 

c  COS0 
w 


(3.38) 


where  e  =  c/n.  Otherwise  the  roots  are  two  complex  conjugate  pairs. 
The  constants  are  determined  using  the  boundary  conditions.  In 
the  general  case  one  has  to  solve  four  simultaneous  linear  equations. 
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b.  Constant  viscosity  and  forcing 


Now  the  free  drift  solution  (see  section  3.2) 


u 


f 


c,  ,cos  9 

W  VJ 


(3.39) 


where  u  =  tane  +  pjhf/cwcose,  is  clearly  a  particular  integral 
of  Eq. (3.34) .  Thus 

4 

u  =  u,  +  Z  C„exp(X  x)  , 

r  4=1  1 

4  al~Xl 

v  =  vf  +  Z  -  Ciex P(Xdx)  • 

4=1  v.  ' 

Constants  C^  and  the  solutions  (3.40)  were  calculated  with  a 
computer.  The  following  parameters  were  held  fixed: 

h  =  1.5  m,  L  =  100  km,  uL  =  0,  VL=0  • 

Inequality  (3.38)  is  now  valid  for  e  ^  2.1.  This  lies  in  a  rea¬ 
listic  range,  so  both  real  and  complex  characteristic  roots  xfc 
are  possible. 

9  —1 

We  chose  the  standard  viscosities  of  e  =  2*10  kgs  and  n  = 
9—1 

10*  kgs  .  These  are  one  order  of  magnitude  lower  than  obtained 
by  Hibler  and  Tucker  (1979)  for  Arctic  summer  conditions,  and 
one  order  of  magnitude  larger  than  used  by  Lepparanta  (1981a)  in 
an  ice  forecasting  model  for  the  Baltic  Sea.  The  ratio  C/n  was 
taken  as  the  same  as  the  compressive  -  shear  strength  ratio  in 
the  plastic  model  of  Hibler  (1979) .  The  forcing  term  was  set 
equal  to  0.2  Nm  2  and  directed  30°  to  the  left  from  the  positive 
y-axis,  i.e.  30°  off-ice  with  ice  on  the  right.  This  choice  gives 
nearly  edge-parallel  solutions. 

The  standard  case  results  in  v-velocity  about  10%  less  than 
free  drift  at  the  ice  edge  (Fig. 19).  Viscosity  also  slightly 
changes  the  direction  from  free  drift  and  hence  u-velocity  is 
very  small.  In  the  standard  case,  e=2  and  the  characteristic  roots 
X ^  are  complex.  Consequently,  there  is  an  Ekman-type  oscillation 
present  in  the  solution.  This  has  a  small  amplitude  but  still  it 


(3. 40. a) 

( 3 . 40 . b) 


can  be  seen  in  u-velocity.  For  comparison,  low  and  high  strength 
cases  were  solved  (Fig. 19) ;  they  were  defined  as  having  one 
order  of  magnitude  lower  and  higher  viscosities  than  the  standard 
values  (their  ratio  e  was  the  same). 

Removing  the  bulk  viscosity  does  not  change  the  standard  case 
much  because  we  consider  basically  an  edge-parallel  shear  flow 
(Fig. 20).  The  situation  would  be  different  if  we  had  an  on-ice 
forcing.  With  no  shear  viscosity  v-velocity  becomes  larger  than 
the  free  drift  value  and  even  increases  towards  the  MIZ-interior 
ice  boundary  (Fig. 20).  This  is  easy  to  understand  when  looking 
at  the  original  equations  (3.30):  with  n  =  0,  v  needs  to  satisfy 
a  balance  similar  to  free  drift  but  with  less  retarding  force 
through  u  because  u  is  affected  by  bulk  viscosity. 


c.  Constant  viscosity  and  variable  forcing 

For  simplicity,  we  neglect  now  the  terms  due  to  the  earth's 
rotation,  i.e.  we  set  8=0  and  f  =  0.  In  addition,  the  forcing 
is  aligned  with  the  ice  edge.  Then  u  =  0  and  v  must  be  solved 
from 


d2v 

n— j  -  c  v  +  F  =  0  . 

dx^  w 


(3.41) 


This  has  a  general  solution 

v  =  C1exp(^x)  +  C2exp(-Xx)  +  vx  ,  (3.42) 

where  and  C2  are  undetermined  constants,  *  =  / cw/ n  and  v^  a 
particular  integral  of  Eq. (3.41).  We  take  now  a  quadratic 
polynomial  for  F, 

F  =  Fq  +  FiX  +  F2x2  (3.43) 


and  then 


V1  = 


•(F 


w 


2F2* 


-1 


+  F2x 


(3.44) 
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Utilizing  the  boundary  conditions  the  solution  becomes 


v  =  v 


1 


.. .  cosh(Xx)  F1  # sinh [ X (L-x) ] 
V1 '  '  cosh (XL)  c~T  cosh (XL) 


(3.45) 


We  fixed  as  earlier  L  =  100  km.  For  F  three  cases  are  considered: 

F  =  Fq (1-x/L) n*  n  =  0,  1  and  2,  Fq  =  0.2  Nm-2. 

The  free-drift  solution  is  shown  in  Fig. 21c.  It  illustrates 

also  the  form  of  the  forcing  term  since  now  free  drift  equals 

c"1  times  the  local  force.  A  nonzero  viscosity  smooths  the  curves 

especially  near  the  ice  edge  where  the  no  stress  condition  must 

be  satisfied  (Figs.  21a-b) .  Let  us  look  at  the  point  where  the 

speed  has  dropped  by  20%  from  its  value  at  the  edge.  With  qua- 

8  -1 

dratic  forcing,  it  is  11  km  for  n  =  0,  20  km  for  n  =  10  kgs 

9  -1 

and  30  km  for  n  =  10  kgs  .  Thus  to  produce  a  sharp  increase  in 
the  velocity  near  the  ice  edge  in  a  constant  viscosity  model,  one 
must  have  a  very  sharp  increase  in  forcing. 


d.  Variable  viscosity  and  constant  forcing 

As  in  c.,  we  neglect  the  earth's  rotation  and  align  the  forcing 
with  the  ice  edge.  Our  equation  for  v  is  thus 


Iflng)  -  c„v  4.  F  =  0  .  (3.46) 

This  can  be  solved  for  some  spesific  forms  of  n  only.  We  take  a 
quadratic  function 

n  =  nL(x/L)2  ,  (3.47) 


where  nL  is  the  viscosity  at  the  MIZ-interior  ice  boundary.  Eq. 
(3.46)  becomes  now  the  Euler  equation 


x2  * 
dxz 


2x  —  - 

dx 


CWL  *  L2F 
— ~ —  v  +  =—  =  0 


(3.48) 


which  has  a  general  solution 


V  =  CL  XU1+  C2  x*2  +  F/cw 


(3.49) 
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(3.50) 


where  and  C2  are  undetermined  constants  and 

ui,2  "  ±  ^cwl2/t»l  +  1/4  ‘  1/2 

(let  ii ^  take  the  plus  sign)  .  Using  the  boundary  conditions 
we  have 

v  -  {1  -  (x/L)ul>  .  (3.51) 

cw 

Three  cases  are  considered:  low,  standard  and  high  strength, 

g  q  i  q  _  i 

defined  by  nT  =  10  ,  10  and  10  kgs  ,  respectively.  The  MIZ 

"  -2 
width  L  is  as  earlier  (100  km)  and  F  =  0.2  Nm 

The  curvature  of  the  velocity  distribution  across  the  MIZ  may 

have  either  sign,  depending  on  the  magnitude  of  the  limiting 

viscosity  nL  (Fig. 22).  Since  the  viscosity  and  its  gradient 

vanish  at  the  ice  edge,  the  ice  velocity  there  equals  the  free 

drift  value.  In  the  case  of  high  strength,  20%  drop  in  the 

velocity  is  reached  at  1.6  km  from  the  ice  edge.  Thus,  sharp 

jet-like  features  in  ice  velocity  can  be  produced  by  variations 

in  the  strength.  Several  kinds  of  viscosity  distributions  are 

of  course  possible.  However,  it  seems  reasonable  to  assume,  as 

here,  that  the  viscosity  decreases  to  zero  (or  a  small  value) 

when  the  ice  edge  is  approached. 


3.4.  General  viscous  plastic  model 

To  examine  how  a  plastic  ice  rheology  modifies  marginal  ice 
zone  dynamics  we  have  carried  out  a  series  of  simulations  using 
the  dynamic  part  of  Hibler's  (1979)  model.  In  addition,  analytic 
solutions  for  the  equilibrium  plastic  adjustment  case  are  con¬ 
structed.  This  section  is  based  on  a  recent  paper  of  Lepparanta 
and  Hibler  (1984b) ;  a  more  complete  description  will  be  given 
in  a  subsequent  paper  of  theirs.  Here  we  discuss  some  of  the 
main  features  of  the  results. 

The  viscous  plastic  rheology  of  the  Hibler  model  gives  linear 
viscous  creep  for  small  strain  rates  and  plastic  flow  otherwise. 


While  other  shapes  could  be  employed,  the  yield  curve  used  here 
is  an  ellipse  in  the  principal  stress  space  and  the  plastic  flow 
is  determined  according  to  the  associated  flow  rule.  The  rheology 
is  defined  by 

a  =  2n e  +  [  ( c-n) ez  -  jP]  I  (3.52) 

P  =  P*h  exp [-C (1-A) ]  ,  (3.53) 

C  =  \  P/max{[e2  +  E^/e 2]1/2 ,  eQ>  ,  n  =  c/e2  .  (3.54) 

Here  n  and  r,  are  the  shear  and  bulk  viscosities,  P  is  hydrostatic 

pressure  term,  and  £j  and  are  the  strain  rate  invariants 

equal  to  the  sum  and  difference  of  the  principal  strain  rates. 

*  * 

P  ,  C,  e  and  eQ  are  rheological  parameters:  P  is  the  strength 

constant,  C  is  the  reduction  constant  for  opening,  eQ  is  the 

maximum  viscous  creep  rate,  and  e  is  the  ratio  of  compressive  to 

shear  strength.  Without  y-derivatives  the  strain  rate  invariants 

are  reduced  to  e.j.  =  3u/3x,  eI3;  =  |3u/3x|.  As  a  result,  the  stress 

state  can  lie  only  in  a  restricted  region  within  or  on  the  yield 

curve  (Fig. 23).  The  momentum  equation  can  now  be  written  as 

|^((C+n>|S  -  |)  *  »ihfv  +  Twx  +  Tax  -  0,  (3. 55. a) 

k  ("li»  -  »i»fu  +  'wy  +  *a*  *  0  '  (3 . 55. b) 

The  standard  rheology  parameters  for  the  present  calculations  are: 

P  =  104  Nm-2,  c  =  20,  e  =  2,  =  2*10-7  s  C  and  e  are  the 

o  ^ 

same  as  in  the  original  model.  The  constant  P  has  been  the  main 
tuning  parameter  of  the  model  and  has  ranged  from  5*10^  in  the 
original  model  to  2.75*10*  Nm  in  a  recent  work  of  Hibler  and 
Walsh  (1982).  Our  creep  limit  eq  is  one  to  two  orders  of  magnitude 
larger  than  usually  applied. 

The  numerical  scheme  and  the  computer  code  of  the  original  model 
(Hibler,  1979,  1980a)  has  been  employed.  The  grid  size  is  4  km  and 
the  time  step  half  an  hour.  Initially,  the  MIZ  is  60  km  wide,  and 
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the  ice  cover  is  compact  and  1.5  m  thick.  The  thickness 
corresponds  to  a  characteristic  MI Z  value  given  in  Wadhams  et  al. 
(1979).  At  time  zero  a  constant  wind  starts  to  blow.  Its  speed 
is  always  10  m/s,  and  in  different  runs  its  direction  varies 
around  the  whole  circle  at  30°  spacing.  In  the  following 
discussion,  the  direction,  noted  by  if>  ,  is  defined  as  the  wind 

cL 

vector  direction  and  is  counted  counterclockwise  from  the  x-axis. 
The  interior  ice  is  assumed  motionless.  The  geostrophic  ocean 
current  has  been  neglected  in  our  simulations,  i.e.  the  ice  is 
purely  wind-driven. 

The  free  drift  case  is  a  useful  reference  solution  for  con¬ 
sidering  the  importance  of  ice  interaction  in  ice  drift.  Here 
for  a  1.5-m  ice  thickness,  the  free  drift  speed  is  15.9  cm/s 
and  the  direction  is  10.6°  to  the  right  from  the  direction  of 
geostrophic  wind.  As  the  thickness  decreases  to  zero,  free 
drift  becomes  aligned  with  the  wind  and  its  speed  increases  to 
16.8  cm/s. 

a.  Initialization  of  the  MIZ  flow 

The  initial  velocity  distribution  was  calculated  by  solving 
the  momentum  equation  only.  Ice  thickness  and  compactness  were 
therefore  equal  to  their  initial  values,  1.5  m  and  1.0,  over  the 
whole  MIZ.  The  results  show  a  plastic  slip  at  the  left-hand 
boundary  and  a  nearly  constant  ice  velocity  in  the  MIZ  (Figs. 
24-25).  In  this  "nearly  constant"  flow  the  ice  creeps  in  a 
linear  viscous  manner  and  the  creep  limit  restricts  velocity 
differences  between  neighboring  grid  cells  to  less  than  0.8  mm/s. 
As  an  exception,  in  the  case  f  =  330°  the  yield  strength  is  not 
overcome  and  the  plastic  slip  is  missing. 

The  results  indicate  that  the  free  drift  solution  is  a  very 
good  approximation  if  it  is  directed  more  than  about  45°  away 
from  the  MIZ.  On  the  other  hand,  stricltly  on-ice  forcing  gives 
practically  no  motion  at  all.  Between  these  two  extremes  there 
are  wide  regions  of  nearly  edge-parallel  flows.  Comparing  with  the 
free  drift  it  is  seen  that  the  ice  interaction  has  a  rectifying 
effect  on  these  flows.  It  is  an  important  feature  and  is  not 
difficult  to  distinguish  from  observations.  The  reason  for  this 
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rectification  is  that  ice  has  compressive  strength.  There  is  a 
slight  asymmetry  between  northward  and  southward  flow  patterns 
due  to  the  Coriolis  and  oceanic  turning  angle  effects. 

The  basic  dimensionless  parameter  in  our  problem  is 


P  h 


(3.  56) 


where  L  is  the  MIZ  width.  First,  plastic  deformation  is  possible 

it 

only  for  y  1  y  »  where  y  depends  on  the  mode  of  deformation  at 
the  left-hand  boundary  and  on  the  yield  ellipse  elongation  e. 
Second,  in  plastic  flow  the  effect  of  ice  interaction  can  be 
thought  of  as  reducing  the  wind  stress  by  a  factor  of  1  -  y/y  . 
The  values  of  y  for  some  spesific  cases  are: 


uniaxial  extension 
pure  shear 

uniaxial  contraction 


General 
i(/l+e"2  -  1) 

i*"1 

\(Jl+e~2  +  1) 


Standard  (e=2) 
0.  06 
0.25 
1.06 


it 

The  present  work  has  y  =  0.62,  and  we  can  see  from  Fig. 24  that 

* 

the  1  -  y/y  dependency  describes  the  solution  qualitatively  well. 


b.  Adjustment  between  mass  and  velocity 

As  the  MIZ  flow  evolves  with  time,  the  motion  toward  the 
interior  ice  must  come  to  a  stop.  Consequently,  the  MIZ  gets 
separated  from  the  interior  ice,  or  a  strictly  edge-parallel  flow 
develops  (Fig. 26).  The  critical  point  seems  to  be  whether  the 
free  drift  solution  is  inward  or  outward. 

When  separated  from  the  interior  ice,  the  motion  of  the  MIZ 
is  approximately  in  the  free  drift  state.  The  velocity  is  then  not 
exactly  a  constant  but  is  slightly  modified  by  the  Coriolis  effect 
and  the  hydrostatic  pressure  term.  The  variations  are  less  than 
2  cm/s,  which  is  about  one-tenth  of  the  average  speed. 

The  edge-parallelity  is  obtained  through  mutual  adjustment 
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between  the  ice  mass  distribution  and  velocity.  On-ice  forcing 
becomes  balanced  with  the  pressure  gradient  within  the  MIZ.  After 
48  hours  model  simulation  edge-parallel  flows  are  already  well 
developed  and  have  approximately  constant  velocities  except  for 
the  last  grid  cell  (Fig. 26).  The  high  edge  velocities  are 
partially  a  numerical  artifact  due  to  the  fixed  grid.  One  should, 
however,  note  that  the  velocity  of  the  MIZ  is  generally  lower 
than  the  free  drift  velocity.  Introducing  perturbations  at  the  ice 
edge  would  give  a  free-drifting  ice  edge  and,  hence,  a  sharp 
increase  in  ice  velocity  close  to  the  edge.  The  form  of  the 
compactness  profile  is  very  sharp  (Fig. 27).  For  along-edge 
winds,  the  compactness  drops  from  0.8  to  0  in  8  km  (two  grid 
cells)  . 

We  can  gain  a  perspective  on  the  numerical  solution  by 
analytically  examining  the  steady  state  adjustment  case.  In 
steady  edge-parallel  flow,  the  u-velocity  must  be  identically 
zero.  The  continuity  equation  is  then  automatically  satisfied 
(see  section  3.1)  and  the  momentum  equation  becomes 

~  \  a£  +  pihfv  +  pwCwsin0'v'v  +  Tax  s  0  '  (3. 57. a) 


dv 

3x 


} 


-  Pwcwcos0[v|v  +  Tay  =  0  ,  (3. 57. b) 


and  now 


n  =  P/{2e2  max(EQ,  e_1|gj|)}  .  (3.58) 

Since  P  =  P(h,A),  our  problem  contains  three  unknowns:  ice 
thickness,  compactness  and  velocity.  We  have  only  two  equations 
(3.57).  The  results,  however,  depend  weakly  on  the  thickness  and, 
hence,  we  take  it  as  constant  in  the  following  considerations. 

Without  shear  strength  ice  velocity  must  be  constant,  with  the 
absolute  value  (  I  I/pwC.,cos9  }  '  .  For  finite  shear  strength, 

dy  w  w 

the  ideal  plastic  case  is  approached  as  eQ  -*■  0.  Then,  in  plastic 
deformation 
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dx 


MT  / 1-*  V  v 

=  le  sgn(d^ 


(3.59) 


and  the  pressure  gradient  can  be  easily  eliminated  from  Eq. 
(3.57).  Then  we  have  a  simple  quadratic  equation  for  v,  and  the 
solution  is 


|V|  =  J  i'-U  (i)2  -  jf- 


where 


lTl  =  I Tay I”  e‘1lTax>  ' 

C*  =  p  c  ,(cos8  +  e  "'"sine  ) 

W  W  — 


r  = 


+  e  1pihf  , 


(3.60) 


(3. 61. a) 
( 3 . 61 . b) 
(3.61.C) 


and  in  +  pairs,  +  is  for  southward  flow  and  -  is  for  northward 
flow.  For  a  given  spatially  constant  wind  stress  the  solution  is 
constant.  Thus  plastic  slip  occurs  at  the  left-hand  boundary,  and 
in  the  MIZ  ice  velocity  is  the  given  constant  while  the  stress 
state  remains  indeterminate. 

Next, since  v  is  constant,  it  follows  directly  from  Eq. (3.57.a) 
that  P  is  linear: 


P  =  b(x  -  Lq)  ,  (3.62) 

where  Lq  is  the  location  of  the  ice  edge  after  adjustment  and 

b  =  rax  +  PwCwsine|v|v  +  p^fv  .  (3.63) 

Naturally,  for  the  solution  to  be  valid  b  must  be  positive. 
Through  Eq. (3.53)  we  can  obtain  the  compactness  profile: 

A  =  1  +  C'1  log{-£-  (x  -  L  )  >  .  (3.64) 

P  h  ° 


This  analytic  form  gives  a  very  sharp  profile  (Fig. 27).  For  example. 


if  b  =  0.05  Nm-^  then  A  =  0.78  at  4  km  from  the  edge. 

The  model  calculations  at  48  hours  agree  rather  well  with  the 
ideal  plastic  steady  state.  In  some  cases  the  model  solution  has 
not  yet  fully  developed  to  the  steady  state.  The  applied  creep 
limit  is  so  small  that  the  results  should  really  approximate 
well  the  ideal  plastic  flow.  On  the  other  hand,  if  the  creep 
limit  is  increased,  linear  viscous  law  is  approached.  Then  the 
viscosity  would  be  spatially  varying  and,  as  shown  in  section  3.3, 
it  would  be  possible  to  obtain  sharp  changes  in  ice  velocity 
close  to  the  ice  edge. 


4.  FIELD  STUDY  ON  MIZ  ICE  KINEMATICS 

4.1  MIZEX-83  experiment  in  the  Greenland  Sea 

The  marginal  ice  zone  in  the  Greenland  Sea  has  been  selected 
as  the  region  of  main  field  effort  of  the  MIZEX-program  (Wadhams 
et  al.,  1981;  Johannessen  et  al . ,  1983a).  The  first  experiment 
took  there  place  in  the  summer  1983  and  was  named  as  "MIZEX 
Greenland  Sea  Summer  Pilot  Study";  here  it  is  briefly  called 
"MIZEX-83."  The  experiment  lasted  about  two  months.  Three 
research  vessels  were  used  and  remote  sensing  studies  were 
carried  through  from  satellites  and  aircraft.  Further  general 
information  of  MIZEX-83  is  given  in  Johannessen  et  al.  (1983a). 

The  present  writers  participated 
in  MIZEX-83  during  June  9  to  July  12  on  the  Norwegian  research 
vessel  Polarbjorn.  During  the  experiment  the  Arctic  sea  ice 
edge  lay  about  zero  meridian  west  of  Svalbard  and  between  the 
latitudes  of  80°N  and  8l°N  north  of  Svalbard  (Fig.  13).  Oceanographic 
instruments  were  first  deployed  in  the  Fram  Strait  and  then 
the  research  area  was  selected  for  the  so-called  drift  phase 
when  our  ice  kinematics  measurements  were  to  be  made.  During 
this  phase  the  ship  was  moored  to  an  ice  floe  thus  drifting 
with  the  ice;  it  began  at  0800Z  on  June  28  and  ended  at  1200Z 
on  July  8.  In  the  beginning  the  floe  was  situated  at  81°  03'N 
6°E,  about  20  km  from  the  ice  edge  (Fig.  13). 

In  the  presentation  that  follows  time  is  generally  counted 
as  a  MIZEX  day  which  has  the  zero  point  at  27  June,  00Z.  Thus 
MIZEX  day  transforms  to  Julian  day  by  adding  178. 

The  ice  conditions  were  almost  ideal  for  our  study.  Ice 
concentration  was  seven  to  nine  tenths  around  the  research 
area  and  down  to  the  ice  edge,  with  some  patchiness  close  to 
the  edge.  Horizontally  the  ice  cover  was  very  uniform  and 
floe  diameters  were  typically  around  50  meters.  In  our  measurement 
scaale,  5  km,  the  pack  well  resembled  a  two-dimensional  granular 
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media  (Fig.  28).  The  ice  floes  were  a  mixture  of  first-year 
and  multiyear  ice  with  mean  thickness  close  to  2  m.  During 
the  drift  phase  Polarbjorn  drifted  about  50  nautical  miles 
northeast  which  is  opposite  to  the  long  term  mean  flow  in  this 
region  (Fig.  29) . 

Our  ice  kinematics  measurements  were  made  using  the  Del 
Norte  microwave  trisponder  system.  It  has  been  used  earlier 
by  CRREL  scientists  off  the  north  coast  of  Alaska  (Tucker  et 
al.,  1980),  and  a  similar  method  has  also  been  used  for  ice 
deformation  studies  in  the  Baltic  Sea  (Lepparanta,  1981b). 

The  system  determines  distances  between  its  stations  by  measuring 
the  time  the  microwave  signal  (frequency  about  9.5  GHz)  needs 
to  travel.  The  absolute  accuracy  is,  according  to  the  manufactor, 
three  meters. 

Our  system  configuration  consisted  of  a  master  station 
(two  master  units),  slave  station  and  four  remote  stations 
(Fig.  30).  The  direct  distances  from  the  master  to  the  remotes 
and  the  range  loops  master-slave-remote-master  were  measured. 

One  remote  station  was  located  at  the  same  site  as  the  slave 
station  in  order  to  get  the  baseline  (master-slave)  length 
which  is  needed  for  triangulation.  As  can  be  seen  from  Fig. 

30,  the  three  measurements  (dQ,  dj,  sj)  for  each  remote  give 
their  positions  in  the  coordinate  system  aligned  with  the  baseline 
and  the  master  station  as  the  origin  (the  master  was  located 
at  Polarbjorn).  That  is,  the  Del  Norte  system  does  not  provide 
information  of  the  location  and  orientation  of  our  measurement 
array  in  the  earth  coordinate  system.  The  latter  was  measured 
less  accurately  using  radar  reflectors  on  the  remote  sites 
and  the  ship's  gyrocompass  but  the  former  is  not  so  important 
to  us  since  we  are  studying  differential  motion  in  the  ice 
pack.  The  location  was,  however,  routinely  determined  in  MIZEX-83 
with  the  ship's  satellite  navigation  system. 

Differential  motion  can  be  derived  from  position  time 
series.  It  was  initially  planned  to  make  Del  Norte  measurements 
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every  hour,  but  once  data  began  to  be  acquired  it  became  obvious 
that  shorter  time  intervals  would  be  necessary  to  resolve  all 
motion  occurring.  Accordingly  the  acquisition  computer  was 
reprogrammed  to  automatically  acquire  data  at  3  min  intervals. 

This  operational  mode  was  employed  throughout  the  experiment. 

The  orientation  was  measured  initially  every  two  hours  but 
approximately  at  half-way  of  the  drift  phase  the  interval  was 
changed  to  one  hour. 

In  terms  of  the  overall  field  operation  there  were  no 
major  problems  and  the  system  operated  much  better  than  expected. 
There  are  a  few  short  gaps  in  the  data  mainly  due  to  delays 
in  changing  batteries  in  the  remote  stations  when  the  weather 
was  so  foggy  that  the  helicopters  could  not  fly.  Battery  lifetime 
were  about  4  days,  but  varied  with  different  types  of  transponders 
Gaps  appeared  also  due  to  polar  bear  problems  and  interference 
from  Motorola  radar  transponders.  After  about  5  days  one  master 
unit  began  malfunctioning  but  fortunately  we  had  a  backup  for 
this  particular  unit  and  it  was  replaced. 

The  stations  on  the  ice  were  mounted  on  the  top  of  3m 
(remotes)  or  4.5  m  (slave)  long  pipes  which  stood  with  the 
help  of  rigid  supports.  These  tripods,  designed  at  CRREL , 
turned  out  to  be  very  good  on  soft  summer  ice.  The  Del  Norte 
system  requires  a  line  of  sight  between  the  stations  where 
the  signal  needs  to  travel.  Under  normal  atmospheric  refraction 
conditions  the  length  of  this  line  in  km  equals  about  4-  ►H  , 
where  H  is  the  sum  of  the  heights  of  the  two  particular  stations 
in  m.  However,  the  shadowing  effect  of  ridges  reduces  the 
effective  tripod  height.  We  had  one  minor  problem  in  that 
the  tripod  height  was  a  bit  low  for  the  initial  slave  site 
C.  To  correct  this  the  slave  unit  was  moved  to  site  D. 

Overall  the  automated  nature  of  the  system  together  with 
its  high  accuracy  proved  most  valuable.  These  features  allowed 
the  acquisition  of  by  far  the  most  detailed  and  temporally 
dense  ice  deformation  measurements  in  existence  to  date.  The 
data  is  presented  in  a  report  of  Hibler  and  Lepparanta  (1984). 


4.2.  General  features  of  MIZ  ice  kinematics 

We  are  concerned  here  with  two  coordinate  system,  Del 
Norte  and  north-east,  which  are  illustrated  in  Fig.  31.  Their 
orientations  differ  by  the  angle  of  8  which  is  a  function 
of  time.  Angles  are  always  given  here  counterclockwise  from 
x-axis.  Since  kinematic  quantities  basically  deal  with  changes 
of  positions  the  accuracy  of  measuring  displacements  (but  not 
absolute  positions)  is  the  most  relevant. 

The  Del  Norte  measurement  triple  (dQ,  dj,  Sj)  for  a  given 
station  (Figs.  30-31)  is  transformed  to  Cartesian  or  polar 
form  through 

xi  =  di  +  (si/d0)[si/2  -  di  -  dQ]  ,  (4.1. a) 

Yi  =  +  /d I  -  x?  ,  (4 . 1 .  b) 

ri  =  dL  ,  (4. 2. a) 

=  arccos  (x^/d^)  .  (4.2.b) 

Note  that  there  are  two  solutions  which  is  clear  for  symmetric 
reasons  (Fig.  31).  Hence  the  sign  of  y-coordinate  must  be 
known.  This  causes  problems  if  a  station  gets  close  to  the 
x-axis.  Then  the  triangulation  error  also  becomes  large  and 
in  fact  using  a  radar  reflector  the  ship's  radar  gives  a  more 
accurate  position  than  the  Del  Norte  system  when  the  direction 
of  the  site  is  within  5  degrees  from  the  x-axis. 

The  displacement  accuracy  of  the  Del  Norte  system  was 
tested  in  the  field  on  several  occasions  by  repeating  the 
measurements  4-5  times  within  a  few  seconds.  In  such  short 
time  the  real  displacement  is  expected  to  be  very  small.  The 
statistics  of  these  data  give  the  relevant  relative  position 
accuracy  and  the  displacement  accuracy  is  then  obtained  by 
multiplying  the  position  accuracy  by  /2  The  accuracy  is 
here  defined  as  twice  the  standard  deviation.  The  result  was 
that  both  the  direct  distance  and  the  range  loop  have  1-m  accuracy. 
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This  value  was  then  used  to  generate  the  accuracy  of  the  Cartesian 
and  polar  coordinates  through  Monte  Carlo  simulation  (Table 
2).  The  site  J  gives  the  worst  readings  because  its  direction 
makes  the  smallest  angle  with  the  x-axis. 

Table  2.  The  position  accuracy  (twice  the  standard  deviation) 
of  the  Del  Norte  Measurements. 


Site 

D 

0 

C 

J 

6x  (m) 

1 

8 

8 

14 

6y  (m) 

- 

10 

3 

22 

6r  (m) 

1 

1 

1 

1 

<5<J>  (deg) 

- 

0.14 

0.11 

0.3 

Radar  reflectors  were  located  on  sites  0,  C  and  J.  For 
them  the  Del  Norte  system  measures  the  angle  $  whereas  the 
radar  gives  the  angle  ?+8  (Fig.  31).  Thus  the  average  difference 
between  the  Del  Norte  and  radar  angles  is  a  good  estimator 
of  6.  Considering  that  the  accuracy  of  the  radar  is  about 
one  degree,  the  averaging  gives  us  the  north-east  orientation 
of  the  Del  Norte  array,  8  ,  to  an  accuracy  of  0.5  deg.  Again, 
Monte  Carlo  simulation  was  applied  to  have  the  accuracies  of 
the  Cartesian  and  polar  coordinates  in  the  north-east  system 
(Table  3). 

Table  3.  The  position  accuracy  (twice  the  standard  deviation) 
in  the  north-east  coordinate  system;  Del  Norte  and  ship's  radar 
data  combined. 


Site 

D 

0 

C 

J 

1>x  (m) 

6 

12 

30 

31 

6y  (m) 

20 

44 

22 

41 

5r  (m) 

1 

1 

1 

1 

5<f>  (deg) 

0.6 

0.6 

0.6 

0.6 
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The  raw  data  is  basically  a  time-series  with  slightly 
varying  interval  (within  a  few  seconds  from  three  minutes) 
and  with  gaps  about  10%  of  the  total  length.  It  was  first 
linearly  interpolated  to  an  equally  spaced  3-minute  time  series. 

Then  a  low-pass  symmetric  filter  with  a  sharp  frequency  cutoff 
designed  by  Hibler  (1972)  was  applied  followed  by  resampling. 
Filtering  necessiates  the  loss  of  n  end  points  from  each  end 
if  the  number  of  weights  is  2n+l.  For  the  filtered  time  series 
below  a  filter  with  241  weights  and  the  transition  band  of 
48  to  60  minutes  has  been  used. 

The  array  configuration  at  approximately  beginning,  middle 
and  end  of  the  drift  phase  is  shown  in  Fig.  32.  A  major  shearing 
and  rotation  occurred  during  the  later  period.  There  was  an 
overall  convergence  through  the  drift  phase  consistent  with 
the  generally  north  eastward  drift  which  would  tend  to  push 
the  ice  edge  back  into  the  main  pack  (Fig.  29).  Note,  however, 
that  there  are  fluctuations  on  top  of  this  trend.  There  is 
some  indication  that  the  shearing  is  greater  than  the  convergence. 

The  direct  distances  had  large  fluctuations  in  the  very 
beginning  and  towards  the  end  of  the  experiment  (Fig.  33). 

The  distance  to  the  site  D  was  about  one  third  of  the  other 
distances  and  thus  it  is  seen  that  the  amplitude  of  the  fluctuations 
scales  approximately  linearly  with  length.  Fig.  34  shows  the 
filtered  area  of  the  triangle  OCJ  at  %  hr.  intervals  over  the 
drift  phase.  The  area  change  shows  a  rapid  opening  and  closing 
early  in  the  experiment  followed  by  a  gradual  convergence  and 
then  a  major  closing-opening-closing  cycle  during  the  last 
several  days  of  the  experiment.  This  last  event  is  accompanied 
by  a  major  shearing  of  the  triangle  (Fig.  32). 

The  triangle  area  in  Fig.  34  is  our  best  time  series  of 
the  changes  in  ice  compactness.  They  are  so  large  over  a  long 
period  that  the  local  deformation  theory  based  on  the  first- 
orderderivatives  of  displacement  does  not  necessarily  yield 
accurate  results.  The  local  theory  is,  however,  good  when 
studying  the  time  rate  of  change  of  the  deformation.  Note 
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that  the  triangle  area  is  independent  of  the  rotation  and  thus 
provided  by  the  Del  Norte  data  only.  We  can  only  present  changes 
in  ice  compactness.  Knowing  the  absolute  value  at  at  least 
one  time  point  would  give  a  reference  value  with  which  the 
scale  on  the  left  in  Fig.  34  could  easily  be  changed  to  the 
compactness  values.  This  should  come  true  when  the  aircraft 
remote  sensing  data  gets  processed.  The  minimum  and  maximum 
triangle  areas  are  about  26  and  36  km^ .  If  the  former  corresponds 
to  compactness  1,  then  the  latter  would  correspond  to  the  compactness 
of  0.72. 

With  regard  to  orientation,  the  baseline  rotated  about 
20  degrees  during  the  drift  phase  (Fig.  35).  This  caused  a 
relative  motion  in  the  array  comparable  to  the  strain  motion 
meaasured  with  the  Del  Norte  system.  At  the  distance  of  5 
km,  20-degree  arc  is  1.75  km  long.  Also  is  shown  in  Fig.  35 
the  rotation  of  Polarbjorn  measured  with  her  gyrocompass. 

The  large  bump  in  the  day  6  is  man-caused  because  of  the  necessity 
to  change  the  ship's  mooring.  The  two  rotation  curves  show 
similar  behaviour  but  the  high  frequency  fluctuations  are  much 
higher  for  Polarbjorn. 

Analysis  of  finer  time  scale  data  shows  there  are  a  large 
variety  of  significant  fluctuations  superimposed  on  the  overall 
trends.  Fig.  36  presents  a  sample  of  raw  measured  distances 
at  3  minute  intervals  showing  fluctuations  over  half  hour  intervals. 
Filtered  distances  are  given  in  Fig.  37  and  large  variations 
on  time  scales  of  a  few  hours  are  seen.  The  fluctuations  in 
these  figures  are  undoubtedly  related  to  floe  bumping  as  well 
as  randomness  in  the  forcing  fields.  However,  it  is  also  possible 
that  they  are  due  to  more  organized  floe  interactions  which 
manifest  themselves  as  kinematic  waves  propagating  across  the 
array  (see  Hibler  et  al . .  1983).  A  simple  example  of  such 
wave  is  the  propagation  of  collisions,  after  pushing  in  one 
end,  through  a  series  of  initially  equally  separated  billiard 
balls . 
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An  interesting  feature  in  the  recordings  was  that  occasionally 
there  was  a  sudden  slip  of  the  order  of  tens  of  meters.  Probably 
due  to  an  overshooting  phenomenon  in  ice  stress  the  distance 
changes  rapidly  and  then  remains  approximately  at  a  constant 
level.  An  example  is  shown  in  Fig.  38  where  a  sample  of  the 
distance  dQ  has  been  plotted.  In  the  days  3  and  4  four  such 
slips  occurred. 
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4.3  Velocity  fluctuations 

In  the  scale  of  the  triangle  the  relative  velocity  fluctuations 
seemed  to  be  typically  of  the  order  of  1  cm/s.  E.  g. ,  for 
a  one-hour  cycle  the  displacement  amplitude  is  only  of  the 
order  of  20m  and  thus  a  very  good  measurement  accuracy  is 
needed  to  study  high  frequency  fluctuations.  An  illustrative 
example  is  given  in  Fig.  39  which  shows  the  stretch  speeds 
(rate  of  change  of  distance)  for  the  lines  ship-site  0  and 
ship-site  J.  They  have  been  obtained  through  first  central 
differening  the  time  series  of  the  interpolated  distances 
and  then  filtering  with  our  standard  low-pass  filter  (transition 
band  48  to  60-minute  periods).  It  is  seen  that  fluctuations 
up  to  8  cm/s  occur.  The  total  variance  of  the  stretch  speed 
was  highest  in  the  site  J  (Table  4) . 

Table  4.  The  variance  of  stretch  speed  (cm^s*^) . 

Site  D  0  C  J 

Variance  0.56  1.95  1.64  5.48 

During  the  days  7  and  8  the  Del  Norte  system  had  a  higher 
measurement  noise  due  to  an  unknown  reason  and  this  period 
cannot  be  used  to  study  frequencies  higher  than  about  one 
cycle  per  hour.  The  data  is  thus  split  into  two  parts  which 
correspond  to  different  flow  regimes  (Fig.  33) :  the  early 
part  (day  3.00  to  day  7.25)  was  relatively  quiet  whereas  large 
variations  occurred  during  the  late  part  (day  8.5  to  day  11.5). 

First  we  shall  study  the  high  frequency  part  with  these  separate 
records  and  then  go  to  lower  frequencies  with  combined  filtered 
data . 

The  method  of  spectral  analysis  is  used  below.  The  program 
used  follows  Hibler  and  LeSchack  (1972).  Briefly,  the  spectral 
density  function  p(A)  of  a  discrete  real  valued  time  series 
is  the  Fourier  transform  of  the  covariance  sequence 
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(4.3) 


p(x)  =  c0  +  r  c  cos(2irXn)  ,  |x|  £  1/2 

n=l 

where  X  is  the  frequency  in  cycles  per  time  interval.  The 
spectral  density  was  here  estimaated  using  the  weighted  covariance 
method  with  Parzen  window  (e.g.  Koopmans ,  1974).  Our  distance 
measurement  error  can  be  considered  as  white  noise  for  which 
p(X)  =  a2,  where  a2  is  the  noise  variance.  For  the  Del  Norte 
system  a2  =(0.5m)2  (Table  1).  It  is  easy  to  see  that  estimating 
the  stretch  speed  u^by  central  differening, 

un(t)  =  [d (t+A t)  -  d(t-At) ] / (2 At)  ,  (4.4) 

leads  to  the  so-called  "blue  noise"  spectrum  for  the  stretch 
speed  measurement  noise, 

Pe  =  1 /2  (a/At)2  {1  -  cos  ( 4 ir X )  }  .  (4.5) 

The  maximum  of  Peoccurs  at  X=l/4  and  then  p  =  (0/At)2 
The  spectra  of  the  interpolated  3-minute  time  series 
is  shown  in  Fig.  40  together  with  the  measurement  noise.  The 
high  frequency  end  falls  off  together  with  the  noise.  The 
signal  equals  noise  at  the  point  where  the  measured  spectrum 
is  10].Og2  ~0.3  above  the  noise  level.  In  all  cases  that  point 
is  at  a  frequency  higher  than  50  cycles  per  day  (about  half-hour 
period).  In  the  frequency  band  from  bout  10  to  50  cycles 
per  day  the  spectral  density  falls  smoothly  with  no  identifiable 
peak  or  change  in  character. 

The  data  was  then  smoothed  with  our  standard  low  pass 
filter  and  the  spectra  was  calculated  for  the  whole  time  series 
record  (Fig.  41).  There  is  some  indication  of  a  peak  near 
the  inertial  frequency  (2.0  cycles  per  day  at  our  location) 
but  the  inertial  motion  is  not  so  clearly  distinguishable 
in  ice  deformation  as  it  is  in  translatory  ice  movement  (cf. 

McPhee,  1978;  Fig.  29  in  this  work).  The  slope  of  the  spectral 
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density  fall-off  seems  to  change  at  about  12  cycles  per  day. 

A  possible  explanation  is  that  at  higher  frequencies  the  randomness 
of  floe-floe  interactions  determines  the  spectra  whereas  at 
lower  frequencies  variations  in  the  external  forcing  fields 
begin  to  dominate  the  variations  of  ice  deformation. 

Small  scale  mechanics  of  pack  ice  should  be  tried  to 
approach  from  a  stochastic  point  of  view.  These  spectra  serve 
as  a  good  verification  data  for  such  studies. 


4.4  Rate  of  deformation  in  the  MIZ 

In  the  length  scale  covering  a  number  of  floes  the  mechanical 
behavior  of  pack  ice  is  commonly  approximated  with  the  frame 
of  a  two-dimensional  continum.  This  approach  is  also  used 
here.  Note  that  the  classical  pack  ice  rheologies  are  based 
on  the  continum  approximation.  Pack  ice  has  no  preferred 
reference  configuration  but  it  has  a  memory  which  enters  the 
rheological  law  through  the  thickness  distribution  function 
(see  E^.  3.3).  In  the  MIZ  short  term  changes  in  the  thickness 
distribution  are  essentially  those  in  the  compactness  which 
was  already  discussed  in  section  4.2.  Thus  this  section  considers 
the  rate  of  deformation,  the  second  independent  variable 
in  the  rheological  law. 

The  rate  of  deformation  is  actually  the  velocity  gradient 
tensor , 
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3v/3x  3v/3y 

This  is  decomposed  into  symmetric  and  anti-symmetric  parts 
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,  (4 . 7.  b) 
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respectively.  The  tensor  e  is  the  strain  rate  tensor  and  w 
is  the  vorticity  tensor  which  has  only  one  independent  component. 
Strain  rate  corresponds  to  physical  deformation,  i.e.  changes 
of  distances  between  the  particles,  whereas  vorticity  corresponds 
to  rigid  body  rotation.  The  former  only  affects  the  ice  rheology 
(see  Eq.  3.3).  Two  useful  invariants  of  the  strain  rate  tensor 
are 

£I  =  trace^)  ' _  (4. 8. a) 

^ll  ^exx  ~  *yy^2+  ^xy  *  (4.8. b) 

The  first  invariant  represents  the  rate  of  change  of  area 
(i.e.  the  velocity  divergence)  and  the  second  twice  the  maximum 
shear  strain  rate. 

It  is  easy  to  understand  that  the  Del  Norte  array,  moving 
with  the  material,  measures  the  physical  deformation,  i.e. 
strain,  exactly.  Thus  the  strain  rate  invariants  can  be  calculated 
using  the  Del  Norte  data  only.  The  individual  components 
of  the  strain  rate  tensor  have  then  no  relevance  since  the 
data  does  not  tell  where  the  coordinate  axes  point  to.  To 
transform  the  strain  rate  tensor  to  the  north-east  system 
and  to  obtain  the  vorticity  one  must  combine  the  Del  Norte 
data  with  the  orientation  observations.  Such  analysis  becomes 
important  when  studying  the  coherence  between  ice  kinematics 
and  the  velocity  fields  in  the  atmosphere  and  the  ocean,  and 
that  will  be  the  next  step  in  our  work  after  the  present  report. 
Note  that  the  rotation  of  the  baseline  (Fig.  35)  does  not 
give  the  vorticity  directly  since  it  is  also  affected  by  shear. 

The  Del  Norte  data  consist  of  very  accurate  distance 
measurements  and  somewhat  less  accurate  angles  (Table  2). 

The  natural  method  to  process  the  strain  rate  tensor  from 
such  data  is  a  certain  kind  of  stretch  distribution  method 
used  much  in  glaciology  and  applied  for  sea  ice  by  Hibler 
et  al.  (1973).  The  method  is  based  on  the  equation  of  the 
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normal  strain  rate  en(i.e.  stretch  speed  divided  with  the 
length  of  the  stretching  line)  in  different  directions  $ 

enU)  =  cos2i})  exx  +  2sin$cos$  exy  +  sin2$  eyy  (4.9) 

Thus  three  normal  strain  rate  measurements  in  known  directions 
give  an  exact  solution  for  the  strain  rate  components.  If 
additional  data  exist,  the  linear  regression  can  be  applied. 

We  used  the  sites  0,  C  and  J  to  obtain  the  strain  rate 
tensor  since  they  correspond  to  the  same  space  scale.  The 
values  of  en  were  obtained  through  central  differencing  the 
interpolated  distances: 

en(t)  =  [a (t+At)  -  d(t-At) ]/[2At-d(t) ]  (4.10) 

and  the  angles  41  were  obtained  as  shown  in  Eq .  (4.2b).  Then 
and  <(>were  smoothed  with  our  standard  low-pass  filter  and 
the  filtered  time  series  was  resampled  every  half  hour.  The 
strain  rate  tensor  was  calculated  from  the  smoothed  record. 

Clearly  the  maximum  shear  rate  fluctuates  with  much  higher 
amplitude  than  the  divergence  (Fig.  42).  In  both  series  there 
is  an  order  of  magnitude  difference  between  the  relatively 
quiet  and  the  more  active  period.  Strain  rates  of  a  few  percent 
per  hour  were  typical  which  is  one  order  of  magnitude  greater 
than  observed  in  the  Central  Arctic  (Hibler  et  al.,  1974) 
or  in  a  compact  ice  cover  in  a  closed  subarctic  basin  (Lepparanta, 
1981b)  . 

Spectral  analysis  gives  a  good  illustration  of  the  difference 
in  shear  and  divergence  (Fig.  43).  The  second  invariant  shows 
6-8  times  larger  variability  than  the  first  invariant  for 
frequencies  higher  than  the  inertial  frequency.  During  the 
active  period  the  inertial  peak  did  not  appear  in  the  maximum 
shear  rate  and  low  frequency  variations  were  very  high. 
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5.  Discussion 

In  this  work  the  dynamics  of  sea  ice  in  the  marginal  ice 
zone  has  been  considered.  Modeling  studies  have  been  carried 
through  with  various  types  of  ice  rheologies  and  field  data 
on  ice  kinematics  obtained  in  MIZEX-83  Greenland  Sea  experi¬ 
ment  are  presented  and  analyzed. 

It  is  recognized  that  the  MIZ  flow  has  certain  specific 
features  under  different  conditions:  1)  Edge  dispersion; 

2)  Banding;  3)  Edge  compaction;  4)  Little  ridging;  5)  Ice 
edge  jet;  6)  Forcing  of  upwelling  in  the  ocean  near  the  ice 
edge.  The  first  and  second  features  occur  in  connection  with 
diverging  ice  edge  and  are  not  much  affected  by  ice  interaction. 
To  understand  and  model  the  rest  of  these  phenomena  a  proper 
formulation  of  ice  rheology  is  required.  It  may  be  possible 
that  the  MIZ  has  viscous  or  plastic  behaviour  depending  on  the 
flow  conditions. 

Ice  edge  dispersion  was  studied  using  the  free  drift  (no 
ice  stress)  approximation.  The  variability  of  ice  thickness, 
such  as  in  Greenland  Sea,  was  shown  to  be  able  to  cause  quite 
notable  differential  drift.  In  particular,  the  thickness  effect 
can  lead  to  floe  clustering  which  was  found  to  be  typical  in 
the  Greenland  Sea  MIZ  in  MIZEX-83.  This  may  be  an  important 
initialization  mechanism  to  the  diffusion  of  MIZ  and  banding. 

Linear  viscous  rheology  was  applied  to  the  MIZ  flow.  It 
is  considered  as  a  rather  crude  approximation  of  ice  rheology 
but,  however,  it  is  much  more  realistic  than  the  free  drift  and 
gives  an  insight  into  how  the  ice  interaction  may  modify  the 
velocity  distribution.  Especially,  it  is  possible  to  obtain 
analytical  solutions  for  the  MIZ  flow  for  certain  interesting 
cases.  With  constant  viscosity  and  forcing,  ice  velocity  does 
not  vary  much  closer  than  30-40  km  to  the  ice  edge.  The  velocity 
at  the  ice  edge  was  107„  lower  than  the  free  drift  for  our 
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standard  viscosity  value.  The  solution  may  contain  a  small 
Ekman-type  spiral.  Monotonical ly  varying  forcing  did  not  have 
a  large  effect  because  the  integral  of  forcing  and  not  the 
local  force  is  most  significant.  (Free  drift  must  follow  the 
local  force  exactly).  On  the  other  hand,  variable  viscosity 
may  produce  quite  different  forms  of  the  velocity  distribution. 
Especially,  a  very  sharp  change  near  the  ice  edge  is  possible. 

As  a  more  advanced  rheology,  Hibler's  (1979)  viscous  plastic 
model  was  applied  to  the  MIZ.  This  model  predicts  a  nearly 
constant  MIZ  flow  under  constant  forcing.  If  the  free  drift 
solution  is  directed  toward  the  open  water  the  MIZ  will  lose 
contact  with  the  interior  ice  and  goes  to  the  free  drift  state. 
Otherwise  there  is  generally  a  plastic  slip  at  the  boundary 
and  the  ice  creeps  elsewhere.  The  ice  velocity  is  naturally 
less  than  the  free  drift  if  the  ice  compactness  is  allowed 
to  equilibrate  in  response  to  a  given  wind  field.  The  profile 
of  ice  compactness  is  sharp  close  to  the  edge.  One  overall 
feature  of  the  model  results  is  that  in  the  absence  of  wind 
nothing  happens,  since  with  a  plastic  rheology  the  ice  will 
not  flow  unless  forced.  This  is  in  contrast  to  Roed  and  O'Brien's 
(1981)  results  where  their  simple  pressure  term  will  cause 
the  ice  to  flow  and  undergo  geostrophic  adjustment  in  the  absence 
of  wind.  Results  also  show  that  divergence  of  ice  near  the 
ice  margin  is  caused  by  ice  interaction.  Hence,  divergence 
of  ice  cannot  be  used  alone  as  an  argument  for  the  absence 
of  ice  interaction. 

Recently  a  considerable  discussion  has  been  devoted  to 
the  so-called  "ice  edge  jet."  Its  existence  has  been  indicated 
by  the  observations  of  Johannessen  et  al.  (1983b)  but  its  nature 
is  not  known.  It  is  clear  that  if  the  MIZ  is  loose  it  has 
no  strength  and  obeys  the  free  drift  law;  the  jet  may  then 
be  caused  by  the  forcing  distribution  only.  For  compact  MIZ, 
the  present  results  show  that  a  jet-like  feature  can  be  present 
as  a  perturbation  in  an  ideal  plastic  model  or  in  a  linear 
viscous  model  with  variable  viscosity. 
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For  future  modeling  work,  two  points  are  recognized  to  need 
a  deep  consideration:  transient  cases  and  stochastic  modeling. 

It  is  not  certain  how  the  velocity  distribution  of  ice  behaves 
during  the  adjustment  process  in  the  viscous  plastic  model  which 
may  take  a  day  or  so.  Probably  this  can  be  solved  using  a 
Lagrangian  grid  only.  We  are  modeling  ice  flow  in  a  rather  small 
scale  where  the  random  component  in  the  velocity  is  large  in 
reality.  Hence,  it  would  be  important  to  know  how  different 
rheologies  react  to  stochastic  forcing. 

It  is  considered  that  more  data  is  needed  on  ice  edge 
dispersion  and  ice  edge  jet,  a  recommendation  to  the  MIZEX-84 
experiment.  In  particular,  it  would  be  relatively  easy  to 
track  the  drift  of  floes  of  different  thicknesses  and  see  if 
the  thickness  effect  is  indeed  significant.  The  nature  of  the 
ice  edge  jet  should  be  shown:  is  it  transient  or  steady,  forced 
or  caused  by  ice  interaction. 

Ice  kinematics  measurements  were  made  with  the  Del  Norte 
microwave  trisponder  system  in  MIZEX-83.  The  study  area 
(about  30km2)  located  northwest  of  Svalbard  about  20 km  from 
the  edge.  Ice  compactness  was  seven  to  nine  tenths  around 
the  area  and  down  to  the  ice  edge.  Horizontally  the  ice  cover 
was  very  uniform.  Floe  diameters  were  typically  arounH  50  m 
and  the  thickness  of  ice  ranged  from  1  to  4  m. 

Our  system  configuration  consisted  of  four  stations  on 
the  ice  around  the  base,  R/V  Polarbjorn  moored  to  an  ice  floe. 

The  Del  Norte  system  measures  positions  of  the  station  relative 
to  the  base  and  thus  can  be  used  to  derive  differential  motion 
in  the  ice  pack.  Our  raw  data  consists  of  the  relative  positions 
at  about  three-minute  intervals.  The  accuracy  of  the  system 
is  very  good.  Distances  are  determined  to  one-meter  accuracy, 
and  the  position  accuracy  is  somewhat  worse  depending  on  the 
triangulation  angles.  The  system  is  the  best  all  weather 
method  available  to  study  small  scale  (order  of  kilometers) 
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ice  kinemetics. 

There  was  an  overall  convergence  of  almost  307,  during 
the  measurement  period.  Towards  the  end  a  major  shear  occurred 
in  the  direction  of  the  ice  edge  and  was  accompanied  by  large 
closing-opening  cycles  in  the  area  of  the  measurement  array. 

On  several  occasions  a  rapid  slip  occurred  in  the  distance 
series  probably  due  to  an  overshooting  phenomenon  in  ice  stress. 

In  the  length  scale  of  about  5  km  ice  velocity  fluctuations 
are  of  the  order  of  1  cm/s.  The  spectra  of  velocity  fluctuations 
showed  a  peak  at  the  intertial  frequency  but  it  was  not  so 
notable  than  generally  in  translatory  ice  drift.  For  frequencies 
higher  than  the  inertial  frequency  the  spectra  had  a  red  noise 
character  down  to  about  half-hour  period  where  the  measurement 
noise  became  nearly  as  large  as  the  signal.  These  were  the 
Tirst  published  high  frequency  spectra  of  ice  motion  and  are 
considered  useful  for  stochastic  modeling  studies. 

The  strain  rate  in  the  MIZ  turned  out  to  be  of  the  order 
of  a  few  per  cent  per  hour  when  the  ice  was  active.  This 
is  about  one  order  of  magnitude  greater  than  is  known  to  be 
in  the  Central  Arctic.  In  frequencies  higher  than  the  inertial 
frequency  the  variance  in  the  maximum  shear  rate  was  about  twice 
as  large  as  that  in  divergence.  At  lower  frequencies  the  shear 
rate  varied  much  more  than  the  divergence. 

The  large  deformation  rates  observed  in  the  ice  pack  are 
probably  of  the  same  order  of  magnitude  than  those  in  the  ocean 
currents.  Thus  it  is  necessary  to  measure  ice  deformation 
very  accurately  in  order  to  study  current  fields  with  instruments 
attached  to  ice  floes. 

It  is  considered  very  valuable  to  carry  through  ice  de¬ 
formation  measurements  with  the  Del  Norte  system  also  in  the 
MIZEX-84  Greenland  Sea  experiment.  The  configuration  used 
in  'TZEX-83  is  applicable,  but  if  two  more  units  are  available 
studies  could  be  made  for  two  different  length  scales.  The 
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field  program  of  MIZEX-84  will  last  six  weeks,  and  the  ob¬ 
servations  will  likely  cover  a  larger  variety  of  processes 
than  in  MIZEX-83.  In  particular,  in  MIZEX-83  the  very  active 
ice  deformation  period  lasted  only  three  days,  and  it  would 
be  useful  to  obtain  more  such  data. 

Overall  the  present  data  set  represents  a  unique  set  of 
measurements  of  the  deformation  field  of  a  series  of  inter¬ 
acting  floes  floating  in  the  ocean.  Because  of  the  fine 
temporal  resolution,  accuracy  of  the  measurements,  and  co¬ 
incidence  of  spatial  current  measurements,  this  data  should 
greatly  aid  in  sorting  out  the  physics  of  marginal  ice  zone 
dynamics  and  kinematics.  Such  efforts  are  currently  in 
progress . 
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Figure  1. 

A  schematic  picture  of  the  marginal  ice  zone  (MIZ)  in  the 
North  Atlantic  on  March  14,  1978.  The  ice  edge  and  fast 
ice  region  are  based  on  U.S.  Navy  ice  charts.  (Weeks,  1980) 


Figure  2. 

A  schematic  diagram  of  the  subzones  in  the  MIZ  (Bauer  and 
Martin,  1980). 


Figure  3. 

The  mean  and  extreme  ice  draft  profiles  across  East  Greenland 
Polar  Front,  October  1976.  (Wadhams  et  al. ,  1979) 


Figure  4. 

Frequencies  of  keels  and  sails  from  transects  across  the 
Denmark  Strait.  (Kozo  and  Diachok,  1973) 


Figure  5. 

Distributions  of  sail  height  and  keel  depth  across  the  Denmark 
Strait.  (Kozo  and  Diachok,  1973) 


Figure  6. 

Ice  drift  in  the  Fram  Strait  based  on  Landsat  images: 

(a)  21  April-7May  1976,  northerly  wind  5-15  m/s  for  most  of 
the  period;  (b)  5-16  May  and  1-12  June  1976,  mostly  calm  weather. 
(Vinje,  1977) 


Figure  7. 

Daily  average  ice  drift  and  surface  current  from  ARGOS  buoys, 
north  of  Svelbard,  1979.  (Johannessen  et  al.,  1983b) 


Figure  8. 


Daily  strain  of  the  ice  cover  in  the  northern  Baltic  Sea  based 
on  LANDSAT  images.  (Lepparanta,  1982) 

Figure  9. 

A  schematic  illustration  of  the  three  classical  pack  ice  rheologies. 

Figure  10. 

A  schematic  illustration  of  the  rheology  of  granular  materials. 

Figure  11. 

Frequency  (in  tenths)  of  sea  ice  concentration  above  4/10  at 
the  end  of  June,  1971-1980.  (Vinje,  1981) 


Figure  12. 

A  contoured  map  of  the  mean  ice  draft  in  Fram  Strait 
22  April  -  4  May,  1979.  (Wadhams,  1983b) 


Figure  13. 

Ice  conditions  in  the  beginning  of  the  drift  phase  of  MIZEX-83 
(Horsk  Meteorologisk  Institutt,  1983).  The  location  of  R/V 
Polarbjorn  is  shown  by  ■. 


Figure  14. 

Visual  observations  of  the  MIZ  ice  during  one  transect  in 
MIZEX-83. 

Figure  15. 

The  idealized  MIZ  for  modeling  studies. 
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Figure  16. 

The  function  F  in  the  free  drift  solution  for  quadratic  water 
stress. 


Figure  17. 

Free  drift  solutions  for  a  quadratic  water  stress:  wind  stress 
is  aligned  with  the  y-axis  and  the  dots  mark  the  tips  of  the  ice 
drift  vectors;  ice  thickness  =  0,  1,  2,  4,  8,  and  16m  and  wind 
speed  =  5  and  10  m/s. 


Figure  18. 

Monte  Carlo  simulation  of  floe  clustering.  In  each  case  the 
lines  follow  each  other,  forming  a  single  long  line.  Each 
digit  gives  the  thickness  of  an  ice  floe  rounded  in  meters. 
The  total  time  is  100  minutes  times  the  floe  size  in  meters. 
A  uniform  wind  of  10  m/s  was  assumed. 


Figure  19. 

Standard,  low  strength  and  high  strength  solutions  for  constant 
viscosity  and  forcing.  The  symbols  on  the  right  show 

free-drift  values. 


Figure  20. 

Solutions  with  no  shear  viscosity  and  with  no  bulk  viscosity 
compared  with  the  standard  case.  Viscosity  and  forcing  are 
constant.  The  symbols  on  the  right  show  free-drift  values. 


Figure  21. 

Constant  viscosity  and  variable  forcing  cases:  (a)  standard 
(n  -  109  kgs-i),  (b)  low  strength  (n  -  10® kgs-1,  and  (c)  free 
drift.  The  symbols  on  the  right  show  free-drift  values 

at  the  ice  edge.  The  free  drift  case  also  illustrates  the  form 
of  the  forcing  term. 


Figure  22. 

Variable  viscosity  cases.  The  symbol 


shows  free-drift  value. 


Figure  23. 


The  yield  ellipse.  Thick  lines  indicate  the  possible  stress 
states  in  the  M1Z  model . 


Figure  24. 

The  initial  MIZ  flow  in  the  viscous  plastic  model.  The  degree- 
numbers  on  the  right  give  the  direction  of  the  geos trophic  wind. 


Figure  25. 

The  stress  state  for  the  initial  MIZ  flow.  The  open  circles 
show  the  stress  at  the  left-hand  boundary  and  the  line  — — 
gives  the  stress  through  the  MIZ  for  300°-wind  case. 


Figure  26. 

The  MIZ  flot  after  48  hours  in  the  viscous  plastic  model.  The 
degree-numbers  in  the  right  give  the  direction  of  the  geostrophic 
wind . 


Figure  27. 

The  compactness  profile  after  48  hours.  The  line  — gives 
the  ideal  analytic  profile  for  along  edge  wind;  the  shape  but 
not  location  is  relevant. 


Figure  28. 

Photograph  of  ice  floes  near  the  center  of  the  Del  Norte  array. 

The  location  of  Polarbjorn  is  noted,  which  for  scale  is  50  m  long. 
This  photograph  was  taken  by  Dr.  Vernon  Squire  (Scott  Polar 
Research  Institute,  Cambridge,  England)  as  a  part  of  ice  photo¬ 
graphy  survey  carried  out  by  helicopters. 


Figure  29. 

The  path  of  Polarbjorn  during  the  drift  phase;  the  numbers  refer 
to  Julian  day  (McPhee,  1983). 


Figure  30. 


Del  Norte  stations  in  MIZEX-83.  The  site  names  and  measured 
quantities  are  indicated. 


Figure  31. 

Del  Norte  (x,  y)  and  north-east  (xn>  yn)  coordinate  systems. 
The  sites  are  also  indicated. 


Figure  32. 

The  array  configuration  at  approximately  beginning,  middle  and 
end  of  the  experiment. 


Figure  33. 

Filtered  (low-pass,  transition  band  48  to  60  minute  periods)  direct 
distances  to  the  Del  Norte  sites. 


Figure  34. 

Filtered  fl.ow  pass,  transition  band  48  to  60  minute  periods)  area 
of  the  triangle  OCJ. 


Figure  35. 

Direction  of  the  baseline  and  the  ship's  heading. 


Figure  36. 

Variation  of  raw  measured  distances  to  Del  Norte  sites  at 
3-minute  intervals  in  early  part  of  experiment.  The  absolute 
distance  values  at  the  beginning  of  the  sample  are  indicated 
on  the  left. 


Figure  37. 


Variation  of  filtered  (low  pass,  transition  band  48  to  60-minute 
periods)  distances  to  Del  Norte  sites  in  the  late  part  of  the 
experiment.  The  absolute  values  at  the  beginning  of  the  sample 
are  indicated  on  the  left. 


Figure  38. 

Raw  measured  distance  to  the  site  D  during  day  2.00  to  day  6.00. 


Figure  39. 

Rate  of  change  of  distance  between  Polarbjorn  and  site  0  (a) 
and  J(b)  during  30  hours  in  the  late  part  of  the  experiment. 
The  time  series  has  been  filtered  with  a  low  pass  filter 
(transition  band  48  to  60  minutes). 


Figure  40. 

The  high  frequency  spectra  of  the  rate  of  change  of  distance 
between  Polarbjorn  and  (a)  site  D,  (b)  site  0,  (c)  site  C, 
and  (d)  site  J.  The  measurement  noise  is  also  shown. 


Figure  41. 

The  spectra  of  the  rate  of  change  of  distance  between  Polarbjorn 
and  sites  D,C,C  and  J. 


Figure  42. 

The  strain  rate  invariants  in  the  Del  Norte  array  estimated  from 
filtered  data  (low  pass,  transition  band  48  to  60-minute  periods) 
a)  divergence,  b)  twice  maximum  shear  rate. 


Figure  43. 

The  spectra  of  the  strain  rate  invariants  in  the  Del  Norte  array. 
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